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Abstract. We study the asymptotics in n for n-dimensional Toeplitz determinants whose symbols 
possess Fisher-Hartwig singularities on a smooth background. We prove the general non-degenerate 
asymptotic behavior as conjectured by Basor and Tracy. We also obtain asymptotics of Hankel 
determinants on a finite interval as well as determinants of Toeplitz+Ifankel type. Our analysis 
is based on a study of the related system of orthogonal polynomials on the unit circle using the 
Riemann-Hilbert approach. 



1. Introduction 

Let f{z) be a complex-valued function integrable over the unit circle with Fourier coefficients 

1 f 

fj = ^J^ f{e'')e-'^'de, j = 0, ±1, ±2, . . . 
We are interested in the n-dimensional Toeplitz determinant with symbol f{z), 

(1-1) DM{z)) = det(/,_fc)",fcio- 

In this paper we consider the asymptotics of Dn{f{z)) as n ^ oo and of the related orthogonal 
polynomials as well as the asymptotics of Hankel, and Toeplitz+Hankel determinants in the case 
when the symbol /(e*^) has a fixed number of Fisher-Hartwig singularities [23, 33], i.e., when /(e*^) 
has the following form on the unit circle C: 



(L2) f(z) = e^(^)zS.-oft- Yl\z- z,\^^^g,^,p^{z)z-''\ z = 6 G [0,27r), 

j=0 

for some m = 0, 1, . . . , where 

(1.3) Zj = e'^^ , J = 0, . . . , m, = Oq < 9i < ■ ■ ■ < Om < ^tt; 

,^ .N.N fe'^^' < argz < 0^ 

(1-4 g, ,3 {z) = g[^ {z) = ■> ' 

^ ^ ^ I e ^^'^^ 9j < avgz < 27r 

(1.5) ^aj > -1/2, (3j £C, j = 0, . . . , m, 

and y(e*^) is a sufficiently smooth function on the unit circle (see below). Here the condition on 
aj insures integrability. Note that a single Fisher-Hartwig singularity at zj consists of a root-type 
singularity 



(1.6) \z — Zj 



2a, 



2 sin ■ 



2 



and a jump gf^.{z). A point zj, j = 1, . . . ,m is included in (1.3) if and only if either aj 7^ or 
Pj / (or both); in contrast, we always fix zq = 1 even if ao = /3o = (note that gj3g{z) = e"*'^^"). 
Observe that for each j = 1, . . . , m, z^^ gp. [z) is continuous at z = 1, and so for each j each "beta" 

singularity produces a jump only at the point Zj. The factors z - are singled out to simplify 
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comparisons with existing hterature. Indeed, (1.2) with the notation b{9) = eY^^''^^ is exactly the 
symbol considered in [23, 4, 5, 6, 7, 8, 9, 12, 13, 14, 20, 21, 41]. We write the symbol, however, in a 
form with z^i=^^^ factored out. The present way of writing f{z) is more natural for our analysis. 
On the unit circle, V{z) is represented by its Fourier expansion: 

(1.7) V{z)= Vkz', Vk = ^l V{e'')e-^''de. 

k=—oo ^ 

The canonical Wiener-Hopf factorization of e^^^^ is 

(1.8) e^(") = 6+(z)e^''6_(z), h+{z) = e^'^=^^>''\ 6_(z) = e^^=— 

In the case that aj = (3j = 0, f{z) = e^^^\ the classical strong limit theorem of Szego (in its 
most general form, see, e.g., [37]) asserts that as ^ oo, 

(1.9) DnU) = exp jnl^o + ^^fc^-^j (1 + o(l)), 

provided V{z) G H'l^ = {V = EZ-oo Vkz" : EZ-o. MVk\' < oo}. 

Fisher and Hartwig [23] were led to single out symbols of type (1.2) based on the solution of a 
variety of specific problems from statistical mechanics, in particular, the solution of the spontaneous 
magnetisation problem for the Ising model. Indeed the square of the magnetisation can be expressed 
as the limit as n — oo of a Toeplitz determinant Dn{f) (which represents a 2-spin correlation 
function at distance n between spins) where the symbol / is a particular example of (1.2) and has 
the following properties depending on whether temperature T is lower, equal or higher than the 
critical temperature T^. 

• for T < Tc, f has no Fisher-Hartwig singularities; 

• for T = Tc, f has one singularity dX zq = 1 with = 0, (3o = —1/2; 

• for T > Tc, / has one singularity at = 1 with ao = 0, /?o = — 1- 

For / of type (1.2), Fisher and Hartwig made a general conjecture in [23] about the asymptotic 
form of Dnif), 

(1.10) Dnif) ^ En" e"-^", n^oo, 

where a = EJLoi'^'j ~ (^j)^ ^^'^ E is & constant depending on /. Considerable effort has been 
expended in the mathematics and physics communities in verifying (1.10). 
Introduce the seminorm: 

(1.11) 111/3111= max |3fi/3,-3f?/?fe|. 

The case when < 1, i.e., when all lie in a single half-closed interval of length 1, namely 

G (g — 1/2,0' + 1/2]) G has now been essentially settled (see however Remark 1.7 below): 
In [41], Widom proved the conjecture when > —1/2, and all (3j = 0. In [4], Basor then verified 
the conjecture when ?R.aj > —1/2, and = 0. In [12], Bottcher and Silbermann established the 
result in the case that \?fiaj\ < 1/2, |3f?/3j| < 1/2. Finally, in [21], Ehrhardt verified the conjecture 
for ^aj > —1/2, < 1. In these papers, the explicit form of E was also established (see [21] 

for a review of these and other related results) . 
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Theorem 1.1. (Ehrhardt [21]j. Let /(e*^) he defined in (1.2), V{z) be C°° on the unit circle, 
'^aj > —1/2, |3fi/3j- — 3fi/3jt| < 1, and aj =b Pj ^ —1, —2, . . . for j,k = 0,1, ... ,m. Then as n ^ oo, 



(1.12) £»„(/) = exp 



3=0 

X n 

0<j<fc<m 



fc=i 

^ G(l + a,+/?,)G(l + a,-/3,) 
J-J- G(l+2a,) + 



X 

i=o 

where G{x) is Barnes' G-function. The double product over j < k is set to 1 ifm = 0. 
Remark 1.2. The branches in (1.12) are determined in the natural way as follows: b±{zj)~°'J^^- 



exp{(-a, ± Pj) Er=i V±kZ^% (^,z7ie-^-)"^-^'=-"'=^^ = exp{i(0fe - Oj - n){ajPk - akPj)}, and the 



remaining branches are principal. 

Remark 1.3. In the case of a single singularity, i.e. when m = or m = 1, ao = /3o = 0, the 
theorem implies that the asymptotics (1.12) hold for 

(1.13) 3fi!a„>-^, PmeC, am±Pm^-h-2,... 

In fact, if there is only one singularity and y = 0, an explicit formula is known [12] for Dn{f) in 
terms of the G-functions. 

Remark 1.4. If all ^pj G (-1/2, 1/2] or all G [-1/2, 1/2), the conditions aj ± Pj 7^ -1, -2, . . . 

are satisfied automatically as Jftoj > —1/2. 

Remark 1.5. Since G{—k) =0, k = 0,1, . . . , the formula (1-12) no longer represents the leading 
asymptotics if aj + Pj or aj — Pj is a negative integer for some j. A similar situation arises in 
Theorem 1.13 below if some representations in M are degenerate. These cases can be approached 
using Lemma 2.4 below, but we do not address them in the paper. 

Remark 1.6. Assume that the function V{z) is analytic. Then the following can be said about the 
remainder term. If all Pj = 0, the error term in (1.12) is of order 0(n~^lnn). If there is only one 
singularity the error term is also 0{n~^ Inn). In the general case, the error term depends on the 
differences Pj — Pk. Our methods allow us to calculate several asymptotic terms rather than just the 
main one presented in (1.12) (and also in (1.27) below). In [16], we show that the expansion (1.12) 
with analytic V{z) is uniform in all aj, Pj for Pj in compact subsets of the strip \^Pj — ^Pk\ < 1, 
for aj in compact subsets of the half-plane ^aj > —1/2, and outside a neighborhood of the sets 
aj lb Pj = —1, —2, .... It will be clear below that given this uniformity. Theorems 1.20, 1.25 also 
hold uniformly in the same sense, while for Theorem 1.13 one should replace Pj with Pj (see below) 
in the condition of uniformity. 

Remark 1.7. Theorem 1.1 as proved by Ehrhardt (and as a consequence, Theorems 1.13, 1.20, 1.25 
that we proved below) hold for C°° functions V{z) on the unit circle. In [16], we extend Theorem 
1.1 to less smooth V{z). Namely, it is sufficient that the condition 



(1.14) J2 \km\<^ 
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holds for some s (and hence for all values in (0, s)) such that 



(1.15) 



s > 



1 — 2 maxj — lo\ 



where ui is defined in (4.63) below so that 2maxj — lv\ < 1. In the present work, we show that 
given Theorem 1.1 with the condition (1.15) on V{z), Theorems 1.20, 1.25 hold for V{z) under a 
similar condition with m replaced by r + 1 and contributions from ao, ar+i appropriately changed, 
while Theorem 1.13 holds under the condition (1.28) of Remark 1.15 below. The uniformity in a-, 
/3-parameters will also hold provided s is taken large enough. 

In [16], we give an independent proof of Theorem 1.1, in the spirit of [19, 28, 31], using a 
connection of -Dn(/) with the system of polynomials orthogonal with weight f{z) (1.2) on the unit 
circle. These polynomials also play a central role in the proofs presented here. 

It follows, in particular, from Theorem 1.1 that all Dfe(/) 7^ 0, k = ko,ko + l..., for some 
sufficiently large /cq if ctj i Pj 7^ —1, —2, .... Then the polynomials 4>k{z) = Xk^^ + " " " > (t^k{z) = 



(1.16) 



of degree k, k 

2lT 



ko,ko + 1, . . . , satisfying 



(t}k{z)z ^f{z)de = Xk Sjk, 



1 

2^70 



M^~'>'f{^)de = Xk^Sjk, 



j = 0,l,...,k, 



exist. It is easy to see that they are given by the following expressions: 



(1.17) 



(l>k{z) 



fc+i 



/oo 

/lO 

1 



/oi 

/ll 

fk-11 

z 



fok 
flk 



fk-lk 
^k 



(1.18) 



where 



4>k{z ^] 



/oo foi 
fio /ll 



ifc-i 



/o 

/lik-1 z ^ 
fkk-l Z~^ 



fkO fkl ■ ■ 

fst = ^ J^^ f{z)z-^'-'Ue, s,t = 0,1,..., k. 



We obviously have 
(1.19) 



Xk 



D 



fc+i 



These polynomials satisfy a Riemann-Hilbert problem. In Section 4, we solve the problem asymp- 
totically for large n in case of the weight given by (1.2) with analytic V{z), thus obtaining the 
large n asymptotics of the orthogonal polynomials. The main new feature of the solution is a con- 
struction of the local parametrix at the points Zj of Fisher-Hartwig singularities. This parametrix 
is given in terms of the confluent hypergeometric function (see Proposition 4.1). A study of the 
asymptotic behavior of the polynomials orthogonal on the unit circle was initiated by Szego [38]. 
Riemann-Hilbert methods developed within the last 20 years allow us to find asymptotics of or- 
thogonal polynomials in all regions of the complex plane (sec [18] and many subsequent works by 
many authors). Such an analysis of the polynomials with an analytic weight on the unit circle was 
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carried out in [34], and for the case of a weight witli aj-singularities but without jumps, in [35]. We 
provide, therefore, a generahzation of these results. Here we present only the following statement 
we will need below for the analysis of determinants. 

Theorem 1.8. Let /(e*^) be defined in (1.2), V{z) he analytic in a neighhorhood of the unit circle, 
and (j)k{z) = Xk^^ + ■ ■ ■ ; 4>k{z) = Xk^^ + • • • he the corresponding polynomials satisfying (1.16). 
Assume that \^fij — 5R/3jt| < 1, aj ± j3j ^ —1, —2, . . . , j. A; = 0, 1, . . . , m. Let 

(1.20) ^ = m«v^2Sft(/3^-A-l) 

Then as n ^ oo. 



S = maxra"' 



(1.21) xl-i = exp 



Zk 



where 

(1.22) Vj = exp \ 

Under the same conditions, 
(1.23) 



r(l + Qj + /3j)r(l + ak- Pk) b+{zj)b-{zk) 



I'k 



T{aj - Pj)T{ak + (3k) b-,{zj)b+{zk) 

+0{5^) + 0{5/n)) , 



'i-i 



-ZTT 



lP=0 




i.n„ r{l + aj + Pj)b+{zj) 



+ 



n 



n 

max 

fe n 



(1.24) (/.„(0) = Xn 



r(i 



a,- 



r(«,+/3,-) 6+(z,) 



n 



max ■ 
it 



n 



Remark 1.9. The error terms here are uniform and differentiable in all aj, Pj for Pj in compact 

subsets of the strip {"^R-Pj — '?R-Pk\ < 1, for aj in compact subsets of the half-plane > —1/2, and 
outside a neighborhood of the sets aj ± Pj = — 1, —2, .... li aj + Pj = or aj — Pj = for some j, 
the corresponding terms in the above formulae vanish. 

Remark 1.10. Note that the terms with n?il^'^~l^3~'^) in (1.21) become larger in absolute value than 
the 1/n term for \\\p\\\ > 1/2. 



Remark 1.11. With changes to the error estimates, this theorem can be generalized to sufficiently 

smooth V{z) using (1.19), a well-known representation for orthogonal polynomials as multiple 
integrals, and similar arguments to those we give in Section 6.2 below. 

Our first task in this paper is to extend the asymptotic formula for Dn{f) to arbitrary Pj G C, 
i.e. for the case when not all 3f?/3j 's lie in a single interval of length less than 1, in other words for 

> 1. 



We know from examples (see, e.g., [12, 10, 21]) that in general, the formula (1.1) breaks down. 
Obviously, the general case can be reduced to 3fi/3j G (g — 1/2, g -|- 1/2] by adding integers to Pj. 
Then, apart from a constant factor, the only change in f{z) is multiplication with z^ , £ G Z. 
However, as we show in Lemma 2.4, the determinants Dn{f{z)) and Dn{z^ f{z)) are simply related. 
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They differ just by a factor which involves Xk^ <t>k{^)^ (f^ki^) for large k (these quantities are given 
by Theorem 1.8), as well as the derivatives of the orthogonal polynomials at 0. The derivatives can 
be calculated similarly to 0fe(O), 0fe(O). Thus it is easy to obtain the general asymptotic formula for 
Dn{f{z)). However, this formula is implicit in the sense that one still needs to separate the main 
asymptotic term from the others: e.g., if the dimension I of F„ in (2.9) is larger than the number of 
the leading-order terms in (1.23), the obvious candidate for the leading order in F„ vanishes (this 
is not the case in the simplest situation given by Theorem 1.18). We resolve this problem below. 

Following [10, 21], we define a FH-representation of a symbol. Namely, for f{z) given by (1.2) 
replace (3j by Pj + rij, nj G Z if Zj is a singularity (i.e., if either (3j / or / or both). The 
integers nj are arbitrary subject to the condition ^JLq '^i — ^- In a slightly different notation 
from [10, 21], we call the resulting function /(z; no, . . . , n^) a FH-representation of f{z). (The 
original f{z) is also a FH-representation corresponding to no = • • • = rim = 0.) Obviously, all 
FH-representations of f{z) differ only by multiplicative constants. We have 

m 

(1.25) f{z) = llz]' xf{z;no,...,nm). 

j=0 

We are interested in the FH-representations (characterized by {nj)jLQ) of / such that X^jLo(3ft/?j + 
rij)"^ is minimal. There is a finite number of such FH-representations and we provide an algorithm 
for finding them explicitly (sec the proof of Lemma 1.12 below). We denote the set of such FH- 
representations by A^. Furthermore, we call a FH-representation degenerate if aj + + rij) or 
aj — {/3j+nj) is a negative integer for some j. We call M. non-degenerate if it contains no degenerate 
FH-representations. The set M can be characterized as follows. For a given (3 = (/3o, ■ ■ ■ ,Pm) let 
us call 

m 

Of, = 0:Pj = (3j+nj,J2nj = O} 

j=o 

the orbit of p. In other words, it is the set of P corresponding to all the FH-representations of /. 
We have 

Lemma 1.12. There exist only the following 2 mutually exclusive possibilities. 

• 3/3 G Oj3 such that \\\P\\\ < 1. Then such P is unique and and it is the unique element of 

• 3p & Ofj such that \ \\P\\ \ = 1. Then there are at least 2 such P's and all of them are obtained 
from each other by a repeated application of the following rule: add 1 to a Pj with the smallest 
real part and subtract 1 from a Pj with the largest. Moreover, M = {P E Op : \ \\P\\ \ = 1}. 

Proof. Suppose that the seminorm \\\P\\\ > 1. Then, writing pi^^ = Ps + I, P^^^ = Pt — 1, and 
P^P = Pj if j 7^ s,t, where Ps is one of the beta-parameters with 3f?/3s = minj3?/3_,, Pt is one of 
the beta-parameters with ^Pt = maxj^Pj, we see that < \\\P\\\, and / corresponding 

to P^^^ is a FH-rcprcsentation. After a finite number, say r, of such transformations we reduce 
an arbitrary set of Pj to the situation for which either < 1 or = 1. Note that 

further transformations do not change the seminorm in the second case, while in the first case the 
seminorm oscillates periodically taking 2 values, and 2 — Thus all the symbols of 

type (1.2) belong to 2 distinct classes: the first, for which |||/?''''^||| < 1, and the second, for which 
|||/3('')||| = 1. For symbols of the first class, M has only one member with beta-parameters P^^'\ 
Indeed, writing bj = ^Pj, if —1/2 < b^J'' — q < 1/2 for some g G M and all j, then for any {kjY-'^Q 
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such that YlT=o = and not all kj are zero, we have 



j=0 i=0 j=0 j=0 j=0 j=0 j=Q 



where the first inequality is strict as at least one kj > 0. For symbols of the second class, we can 
find g G M such that —1/2 < b^p — g < 1/2 for all j. Equation (1.26) in this case holds with 
">" sign replaced by ">". Clearly, there are several FH-representations in M in this case (they 
correspond to the equalities in (1.26)) and adding 1 to one of /ji'^^ with 6^^^ = mmj = q — 1/2 

while subtracting 1 from one of P^^^ with 6^^^^ = maxj b^ = q + 1/2 provides the way to find all of 
them. 

A simple explicit sufficient, but obviously not necessary, condition for Ai to have only one 
member is that all ^l3j mod 1 be different. □ 

In Section 6, we prove 

Theorem 1.13. Let f{z) be given in (1.2), ^aj > -1/2, f5j G C, j = 0, 1, . . . , m. Let M be 
non-degenerate. Then, as n ^ oo. 



(1.27) DM = E n nfi^; no,..., nm)){l + o{l)), 




where the sum is over all FH-representations in M.. Each TZ{f{z; no, • • • , n-m)) stands for the right- 
hand side of the formula (1.12), without the error term, corresponding to f{z;no, . . . ,n^). 

Remark 1.14. This theorem was conjectured by Basor and Tracy [10] based on an explicit example: 
sec Example 1.19 below. The case when the FH-rcprcscntation minimizing EjLo(^'^i 
unique, i.e. there is only one term in the sum (1.27), was proved by Ehrhardt [21]. Note that this 
case is exactly the first possibility of Lemma 1.12. Thus, Theorem 1.13 in this case follows from 
Theorem 1.1 applied to this FH-representation. 

Remark 1.15. This theorem relies on Theorem 1.12 and therefore requires V{z) to be on the 
unit circle. As remarked above, we prove in [16] that Theorem 1.1 holds in fact under the condition 
(1.15). It then follows (see Section 6) that, if has several members. Theorem 1.13 holds for any 

1 + E7=o k^^j? + max |(Ji^,)2, (Ji/3f )2| 
(1.28) s > L 



1 — 2 maxj |3f?/3j — uj\ 

where /3j are obtained from /3j ' by subtracting 1 from all /3j ' with the maximal real part and 
leaving the rest unchanged. The number a; is given by (4.63) below with /3j replaced by 

Remark 1.16. The situation when all aj ± /3j are nonnegative integers, which was considered by 
Bottcher and Silbermann in [13], is a particular case of the above theorem. 

Remark 1.17. The case when all the FH-representations of / are degenerate (not only those in M) 

was considered by Ehrhardt [21] who found that in this case Dn(f) = ©(e^^^n^), where r is any 
real number. We can reproduce this result by our methods but do not present it here. 

We will now discuss a simple particular case of Theorem 1.13 and present a direct independent 
proof in this case. 
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Theorem 1.18 (A particular case of Theorem 1.13). Let the symbol f^{z) be obtained from f{z) 
( 1.2) by replacing one fSj^ with /Sj^zizl for some fixed < jo < m. Let^aj > -i, 3?/3j G (-1/2, 1/2], 
J = 0, 1, . . . , m. Then 

(1.29) DM^iz)) = zr-^DMz)), i^n(r(^)) = zl^DMiz)). 

Xn Xn 

These formulae together with (1.23,1.24,1-21,1.12) yield the following asymptotic description of 
Dn{f^)- Let there be more than one singular points zj and all aj it /3j ^ 0. For f'^{z), let (3j^, 
p = 1, . . . ,s be such that they have the same real part which is strictly less than the real parts of 
all the other j3j, i.e. ^Pj^ = ■ ■ ■ = 'Rf3j^ < mmjyj^_...j^ 3fi/3j. For f~{z) let one Pj^, p = 1, . . . , s be 
such that ^Pj^ = • • • = ^/3j^ > maxj^j^^...j^ 3ft/3j. Then the asymptotics of Dn{f^) are given by the 
following: 

(1.30) D„(/+) = zT- ^ z^7^,,,+ (l + o(l)), DM-) = zli2 ^r^>.-(l + '^(l))' 

p=i p=i 

where TZj,± is the right-hand side of (1.12) (without the error term) in which Pj is replaced by 
Pj ±1, respectively. 

Proof. For simplicity, we present the proof only for V{z) analytic in a neighborhood of the unit 
circle. Consider the case of f~{z). It corresponds to one of the Pj shifted inside the interval 
(-3/2,-1/2]. Since 

we see that 

f-[z) = -ZjoZ-^f{z). 
Therefore, using the identity (2.12) below, we obtain 

Dniriz)) = {-ZjXDn{z-'f{z)) = z^^DMiz)). 

Xn 

If, for some ji, j2, . . . , jg, we have that ?R.Pj-^ = . . . = ^Pj^ > maxj^jj^,, ^^ ^Pj, then we see from 
(1.24) that only the addends with n^^-'i^^, . . . , n'^f^is-^ give contributions to the main asymptotic 
term of Dn{f~{z)). Using Theorem 1.1 for Dn{f{z)) and the relation G(l + x) = T{x)G{x), we 
obtain the formula (1.30) for Dn[f~{z)). The case of f^{z) is similar. □ 

Example 1.19. In [10] Basor and Tracy considered a simple example of a symbol of type (1.2) for 
which the asymptotics of the determinant can be computed directly, but are very different from 
(1.12). Up to a constant, the symbol is 

(1.31) /(-)(e'^) = |:^' . 

I i, vr < y < zvr 

We can represent /'^'^"^•* as a symbol with /^-singularities Pq = 1/2, Pi = —1/2 at the points zq = 1 
and zi = — 1, respectively: 

(1.32) /(^^)(z) = 51,1/2 (^)5-i,-i/2(^)e^'^/' 

We see that f^^^\z) = f~{z) and jo = 1- Therefore by the first part of Theorem 1.18, we have 

DM^^^\z)) = {-ir^DM{z)l 

Xn 
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where ct>n{z), Xn, Dn{f{z)) Correspond to f{z) given by (1.2) with m = 1, zq = 1, zi = e*"^, 
Po = Pi = 1/2, ao = ai = 0. 

Observing that s = 2, ji = jo = 1 and j2 = and using (1.30) we obtain 

Dnif^'^'^H^)) = (-l)"((-l)"7^l,_ + 7^o,-). 

Since 7^l _ = 7^o - = i2n)-'^/^G{l/2fG{3/2f{l + o(l)), we obtain 

(1.33) Dn{f^^^\z)) = l±^y|G(l/2)2G(3/2)2(l + o(l)), 

which is the answer found in [10]. 

As noted by Basor and Tracy, f^^^\z) has a different FH-representation of type (1.2), namely, 
with (3q = —1/2, /?! = 1/2, and we can write 

(1-34) f^^^\z) = -9,,-i/2iz)g-i,i/2iz)e-'^/'. 

This fact was the origin of their conjecture. In the notation of Theorem 1.13, the symbol (1.32) 
has the two FH-representations minimizing Yl]=oi^^j + %)^' with no = ni = and the other 
with no = —1, ni = 1. 

Note that in the case = we can always assume that 3?/?^ G [—1/2, 1/2]. The beta- 

singularities then are just piece-wise constant (step-like) functions. This case is relevant for our 
next result, which is on Hankcl determinants. 

Let w{x) be an integrable complex-valued function on the interval [—1,1]. Then the Hankel 
determinant with symbol w{x) is given by 



(1.35) Dn{w{x)) = det (^j x^+^w{x)dx^ 

Define w{x) for a fixed r = 0, 1, . . . as follows: 



n-l 



r+1 

(1.36) w{x) = e^(^) n 1^ ~ ^jf "'^j(^) 

l = Ao>Ai>--->A,+i = -l, a;,(x) = r_.^^. '^^^t^'^^^, G (-1/2, 1/2], 



^l^i Rx < Xj 

1 



Po = Pr+1 = 0, ^aj>--, j = 0, 1, . . . , r + 1. 

where U{x) is a sufficiently smooth function on the interval [—1, 1]. Note that we set /?o = Pr+i = 0, 
5R/3j G (—1/2, 1/2] without loss of generality as the functions uJo{x), uJr+i{x) are just constants on 
(—1, 1), and LOj{x; Pj + kj) = {—l)''^uj{x; Pj), kj G Z. 
In Section 7, we prove 
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Theorem 1.20. Let w{x) be defined as in (1.36) with G (—5,5); j = 1,2, ...,r. Then as 
n — 00, 

(1.37) Dniw) = £)^(i)e[('*+"o+"'-+i)^o-"o^(i)-"'-+i^(-i)+5 Er=i '^^fc] 

r 

Y\ b+ {zj ) b- {zj ) ""-^ ^'^^ X et^^^""*""^^ ^^=1 "^''^ Eo<i<fc<r+i ("i/^fe )] 



r 

X 



X 

0<i<fe<r+l 



n 1^.- - Afe|-2("^"'=+/^^/5'=) |a,A, - 1 + ^(1-A2)(1-A2) 

;fe<r+l 

1 TTn _ ^2^-(a?+/3^)/2 G'(l + + /3j)G(l + aj - Pj) 



G(l + 2ao)G(l + 2a^+i) G(l + 2^^ 

where A = Ylk=o '^k, V{e^^) = U (cos 0), Zj = e*^-) , Xj = cos ^j, j = 0, . . . , r + 1, and t/ie functions 
b±{z) are defined in (1.8). 

Remark 1.21. -D„(l) is an explicitly computable determinant related to the Legendre polynomials 
(it can also be written as a Selberg integral) , c.f . [42] , 

iim Dm -2-^17 _ 7r"+V^G(l/2)^ 

(1.38) 15n(l)-2 JJ[^ _ ^-^-^^-j^ (1 + o(l)) . 

To prove Theorem 1.20 we use the fact that w{x) can be generated by a particular class of 
functions f{z) given by (1-2). Namely, we can find an even function f of 9 (/(e*^) = f{e~^^), 
9 e [0,2n)) such that 

(1.39) ^ = cos^, a;G[-l,l]. 

We must have (see Section 7 below) that m = 2r+l, 9o = 0, 9r+i = vr, 9m+i-j = 2TT — 9j, j = 1, . . . r. 
If we denote the beta-parameters of f{z) by l3j, we obtain /?o = Pr+i = 0, /?j = —/Sm+i-j = —l3j, 
j = 1, . . . ,r. In particular, Yl]Lo (^3 = ^ ^ remarked above. 

In Section 7 we obtain Theorem 1.20 from Theorem 1.1 and the asymptotics for the orthogonal 
polynomials on the unit circle with weight f{z) using the following connection between Hankel and 
Toeplitz determinants established by Theorem 2.6 below: 

(1.40) D^in.{x)) = 
where w{x) and f{z) are related by (1.39). 

Remark 1.22. Asymptotics of a Hankel determinant when some (or all) of fij have the real part 

1/2 can be easily obtained. For the corresponding f{z) this implies that certain lkf3j = — 1/2 and 
"^Pm+i-j = 1/2 and the rest 5ft/9fe G (—1/2,1/2). Thus, Theorem 1.13 can be used to estimate 
D2n{f{z)). For the asymptotics of 02n(^) in this case we need an additional "correction" Ri term 
(given by (4.69) below) which is now 0{n~'^^^~^) = 0(1). 

Remark 1.23. One can obtain the asymptotics of the polynomials orthogonal on the interval [—1, 1] 
with weight (1.36) by using our results for the polynomials 4>kiz) orthogonal with the corresponding 
even weight on the unit circle and a Szego relation (Lemma 2.5 below) which maps the latter 
polynomials to the former ones. 
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Remark 1.24. Asymptotics for a subset of symbols (1.36) which satisfy a symmetry eondition and 
have a certain behaviour at the end-points ±1 were found by Basor and Ehrhardt in [6]. They use 
relations between Hankel and Toeplitz determinants which are less general than (1.40) but do not 
involve polynomials. For some other related results, see [26, 32]. 

Our final task is to present asymptotics for the so-called Toeplitz+Hankel determinants. We 
consider the four most important ones appearing in the theory of classical groups and its applications 
to random matrices and statistical mechanics (sec, e.g., [2, 25, 30]) defined in terms of the Fourier 
coefficients of an even / (evenness implies the matrices are symmetric) as follows: 

(1.41) dct(/,_fc + /,+fc)--io, det(/,_, - /,+,+2);,fcio' detifj.k±fj+k+i)]jio- 

There are simple relations [39, 29, 2] between the determinants (1.41) and Hankel determinants 
on [—1, 1] with added singularities at the end-points. These are summarized in Lemma 2.7 below. 
It is easily seen that if f{z) is an (even) function of type (1.2) then the corresponding symbols of 
Hankel determinants belong to the class (1.36). Thus a straightforward combination of Lemma 2.7 
and Theorem 1.20 (aided by formulae of Section 7) gives the following 

Theorem 1.25. Let f{z) be defined in (1.2) with the condition /(e*^) = /(e~*^). Let 9r+i = ti" o-nd 
^Pj G (— |, I), j = 1, 2, . . . , r, /3o = Pr+i = 0. Then as n —> oo, 

(1.42) 1)'^+^ = e"^o+5[(°o+"^'-+i+*+*)^o-("="o+*)^(i)-K+i+*)^(-^+Sfcii ^^/?] 

r 

X Jj6+(zj)-"^+^^6_(zj)-"^-^J X e"*''[{""+*+^j"=i"^}^j"=i^^+^i<^<'^<'-(°^'^^~"^^^')] 

l<j<k<r 

X fl zf^'\l - z||-H+/^')|l - 2^.|-2".("0+.)|l + ^.^-2a,{ar+,+t) 

^^K+^^+i+^+W)g(i/2)2 - G(l + a,+/j,)G(l + «,-/3,) 
G(l + ao + s)G{l + ar+i + 1) f}^ G{1 + 2a,) ^ ))' 

where A = ^(ao + ctr+i + s + i) + X^j=i ^''^^ 

(1.43) D'^+'' = detifj^k + fj+k)];^Lo, with p = -2n + 2, s = t = ~ 



(1.44) D^+'' = det{fj.k-fj+k+2)';iLo, with p = 0, s = i = J 



2 

1 . 1 



(1.45) D™ = det{fj-k±fj+k+i)];;lo, with p = -n, s = t^, t = ±^ 
Remark 1.26. For the case ^f3j = 1/2 see Remark 1.22 above. 



Remark 1.27. For the determinant det(/j_jfc + fj+k+i)j^=o in the case when the symbol has no 
a singularities at z = ±1 and < 1/2, the asymptotics were obtained in [7] (see also [8] if 

/ is non-even, aj = 0). Note that for symbols without singularities, i.e. for f(z) = e^^^\ the 
asymptotics of all the above Toeplitz-|-Hankel determinants (and related more general ones) were 
found recently in [9]. 
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The results in this paper make it possible to justify various asymptotic formulae that were ob- 
tained previously in the literature on the basis of the Basor- Tracy conjecture and conjectures on 
Hankel and Toeplitz+Hankel determinants. One example is the probability Psin) of a ferromag- 
netic string of length n in the antiferromagnetic ground state in the XY spin chain. For a certain 
range of parameters, Pe{t^) is given by a Toeplitz determinant with 2 /^-singularities such that 
= 1. Thus, Theorem 1.18 verifies the result on Psin), n — > oo, presented in [22], which the 
authors based on the Basor-Tracy conjecture. Another application arises in the framework of the 
random matrix method in the theory of L-functions. Indeed, the consideration of the mean-values 
of L-functions in orthogonal and symplectic families leads to certain Toeplitz+Hankel determinants 
with a single a-singularity at = 1. Then the application of Theorem 1.25 yields the asymptotic 
formulae which can be used to analyze number-theoretical conjectures in the orthogonal and sym- 
plectic cases [15] in the same way as the Toeplitz formulae of Theorem 1.1 are used in the unitary 
case in [27]. In a similar vein, our results can be used to justify the asymptotic results for correlators 
arising in the theory of the impenetrable Bose gas, that were obtained in [36] on the basis of the 
Basor-Tracy conjecture. 



2. Orthogonal polynomials on the unit circle. Toeplitz and Hankel determinants. 

Here we present aspects of the theory of orthogonal polynomials on the unit circle we use in 
this work. Some of the properties we describe here are well-known (see, e.g. [38], [37]), the others 
not so. We also adapt the theory to complex weights we need in this work, while in the literature 
usually only positive weights are considered. 

Let f{z) be a complex-valued function integrable over the unit circle, and let (pkiz) = XkZ^ + ' ' ' > 
4>k{z) = XkZ^ + • • • , /c = 0, 1, . . . be a system of polynomials in z of degree k with the same for 4>k{z) 
and (pkiz) leading coefficients Xk- These polynomials are called orthonormal on the unit circle with 
weight f{z) if they satisfy (1.16). If f{z) is positive on the unit circle, it is a classical fact that 
Dn{f) 7^ for all n > 1, and such a system of polynomials exists. In general, suppose that all the 
Toeplitz determinants Dn, n = 1,2,... (1.1) are nonzero, Dq = 1. Then the polynomials (pki^) 
and (pkiz) for A; = 0, 1, . . . are given by the explicit formulae (1.17), (1.18) for all k = 1,2, ... . For 
k = set 

(2.1) M^) = Mz) =Xo = ll\fD[. 
Relations (1.16) are then equivalent to 

(2.2) i-^ <^k{zM>^{z-^)j{z)dd = hm. z = e'\ k,m = 0,l,.... 

Thus we constructed the system of orthogonal polynomials under condition that all the Toeplitz 
determinants are nonzero. If we only know that -Dn(/) 7^ for all n > Nq with some Nq > 0, 

then we have the existence of (f)k{z) = XkZ^ H , (t>k{z) = Xkz'' H ) Xit 7^ Oj satisfying (1.16) for 

k = Nq,No + 1,.... 

Remark 2.1. From (1.17,1.18,2.1) we easily conclude: 

a) If f{z) is real on the unit circle, we have ^„(z~-^) = (f)n{z), n = 0, 1, . . . , on the unit circle. 

b) If f{e'') = f{e-''), then M^-^) = M^'')- 
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Lemma 2.2 (Recurrence relations). Let Dn{f) ^ 0, n > 0. The orthogonal polynomials satisfy 
the following relations for n = 0,1, . . . : 

(2.3) XnZ(l)n{z) = Xn+10n+l(2;) - ^n+l(0)2;""^Vn+l(2;~^); 

(2.4) XnZ~^^n{z~^) = Xn+l^n+l{z~'^) - ^n+l(0)^;""~Vn+l(^;); 

(2.5) Xn+lZ'^^niz^^) = Xn$n+liz~'^) - $n+l{0) z'"- (l)n{z) . 

Moreover, 

(2.6) x'+l-Xn = '/'n+l(0)^n+l(0). 

Proof. To prove (2.3) consider the function 

9{z) = Xn<t>n{z) - Xn+lZ~^<i>n+l{z) + 0„+i (O)2;"0„+i (z~^). 

We see that it has zero coefficient at and so g{z) is a polynomial in z of degree n. Therefore 
we can write 

n 

f(^) = ^Ck^k{z), 

where Ck = ^ [q^ 9iz)(l>k{z~^)f{z)d6. This integral is easy to calculate using the orthogonality in 

the form of (1.16) (for example, ^ J^'^ <pn+i{z)z~^4>k{z~^)f{z)d9 = {Xn/Xn+i)6nk), and we obtain 
that all Cfe = 0. Thus g{z) = and (2.3) is proved. 

Similarly, considering gi{z) = Xn(pn{z'^) - Xn+iZ(^n+iiz~^) + 0„+i(O)2;~"(?!)„+i(2;) we show that 
gi{z) = 0, which proves equation (2.4). 

Collecting the coefficients at z"^^ in (2.3) we obtain (2.6). 

Finally, multiplying (2.3) by z~'^~^(f)n+i{0), and (2.4) by Xn+ij adding the resulting equations 
together and using (2.6), we obtain (2.5). □ 

Lemma 2.3 (Christoffel-Darboux identity). Let D^if) / 0, n > 0. For any z, a ^ 0, n = 1,2, . . . , 

n-l 

(2.7) (1 - a-^z) Ma-')Mz) = a-^Ma)z''Mz~') - Ma-')Mz)- 

k=0 

For any z ^ 0, n = 1,2, . . . , 

"-1 ^ ^ / J d ^ \ 

(2.8) 'Y(j)kiz~^)(pkiz) = -n(pniz)(pniz~^) + z i(pniz~^)-^4>niz) - ^„(2:) — 0„(2;"^) j . 

Proof. Consider (1 — a^^z)(6fc(a^^)(/)fc(z), for a fixed /c > 0. Using the recurrence relation (2.3) with 
n = /c to express zcpkiz) in terms of (j)k+i{z) and (j)k+i{z^^), and using (2.5) with A: = n to express 
(/)k{a''^) , we obtain: 

(1 - a~^z)(f>kia~^)(j)k{z) = 4>k{a~^)(t)k{z) - 4)k+i{a~^)(f)k+i{z) 

+ '^^^z^+%^,{z-')a-^-\'^+^^k+,{a-') + ^^a-V.(a)^'+^^-'=-^0fe+i(^) 

Xk+l Xk 

'/'fe+i(o)<?fe+i(o) . („\.,k+i2 (.-^\ 

-a (pk\a)z (pk+i\z ). 



XkXk+l 
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Now expressing in the third summand a^'^^(f)k+i{a ^) from (2.3) with n = k and z = a, and in the 
fourth summand z~''~^(j)k+i{z) from (2.4), and by using (2.6), we obtain 

(1 - a~^z)4)k{a~^)(j)k{z) = 4)k{a~^)(j)k{z) - (})k+i{a~^)(t>k+i{z) 



+ a 



-k-l 



4>k+i{a)z''+^A+i{z ^)-a ^(t>k{a)z^^k{z ^). 



Summing this over k from = to n — 1 yields (2.7). Taking the hmit a ^ z in (2.7) gives 
(2.8). □ 

The next lemma allows us to represent the Toeplitz determinant with symbol z^f{z), where £ is 
any integer, in terms of the one with symbol f{z). 

Lemma 2.4. Let the Toeplitz determinants Dn{f) with symbol f{z) be nonzero for all n> Nq with 
a fixed No > 0. Let $fe(z) = 4>k{z)/xk, ^k{z) = 4>kiz)/Xk, k = Nq, Nq + 1, . . . be the system of 
monic polynomials orthogonal on the unit circle with the weight f{z). Fix an integer £ > 0. Then 
if 



Fk = 

we have 
(2.9) 



^fc(O) 



^fc+i(0) 
■^fc+i(0) 



dz 



i^^kio) £^^k+i{o) 



Dn{z'f{z)) 



^k+i-M 
■ i^^k+e-i{0) 

(-l)'"i^„ 



k = No,No + l,...,n-l, 



-D„(/(2)), H > No. 



In particular, for £ = 1, if 4>k{0) 0, k = Nq, Nq + 1, . . . ,n — 1, we have 

^n</'n(0) 



(2.10) 

Furthermore, if 



Dn{zf{z)) = {-ly 



Xn 



-Dnifiz)), n > No. 



we have 
(2.11) 



$fc(0) $fc+i(0) 

i$k{o) i^k+m ■ 
£^^k{o) £^^k+m ■ 

Dn{z-'f{z)) 



^k+e-M 

£^^fe+^-i(0) 



/O, k = No,No + l,...,n-l, 



-Dnifiz)), 



n > No. 



In particular, for £ = 1, if (pki^) / 0; ^ = -^O; Nq + 1, . . . ,n — 1, we have 



(2.12) 



Dr,iz-'f{z)) = {-l)-^!^DM{z)), n > No. 



Xn 



Proof. We give the proof for £ = 1; the generalization is a simple exercise. Recall that since Dn ^ 0, 
n = A^O) -^0 + 1) • • • ) the polynomials (pn{z) = XnZ^ + . . . , Xn / 0, exist for n = iVo, A'^q + 1, . . . . 

Assume first that Nq = 0. Given the polynomials (pkiz) related to the weight f{z), we will need 
the ones corresponding to the weight zf{z). An analogous construction for polynomials orthogonal 
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(2.13) 



on the real line is known as ChristrofFel's formula (see [38], p. 333). Namely, define qk{z) by the 
expression: 

(pk{z) (pk+iiz) 

<Pk{o) 0fc+i(o) ■ 

We see immediately that qk{z) is a polynomial, and if 0^(0) 0) it has degree k with leading 
coefficient —Xk+i4'k{0)- Moreover, by orthogonality, 

r.27r 

qk{z)z-^zf{z)de = 0, j = 0, 1, . . . , - 1. 



/ 

^0 



27r 



For j = k, 

MO) 0fe+i(o) 

Therefore, for monic polynomials Qk{z) = qk{z) / {—Xk+i<t>k{0)) ■, 

</'fc+i(0) 1 



^ qk{z)z-'zf{z)de 



(t>k+i{0) 



1- Q^{z)z-''zf{z)d9 = -- 



Xk 



<Pk{0) XkXk+i 

Thus, the Toeplitz determinant with symbol zf{z), is given by the expression 



n-l 



(2.14) Dnizfiz)) = l[hk = 



i-iy 



k=0 



<Ao(0) XoXi ■ --Xn-lXr. 



(-1) 



Xn 



Dnifiz)), 



which is equation (2.10). The case of z ^f{z), i.e. equation (2.12), is obtained similarly by 
considering (^^(z"^) instead of (j)k{z). Namely, wc start with the definition 



z ^qk{z ^) = 



Mz ^) ^fe+i(^ ^) 
^fe(O) ^fe+i(0) 

and proceed as before. 

Suppose now that Dn{f) ^ for n > 0, but Fk, Fk are known to be nonzero only for k = 
Nq,Nq + 1, . . . ,n — 1, with some Nq > 0. Consider the polynomials Fk{z), Fk{z) defined as F^, 
Fk with the argument of the orthogonal polynomials replaced with z. Obviously, the set Vt of 
possible zeros of Fk{z), Fk{z) for A; = 0, 1, . . . , iVo — 1 is finite. We now replace (2.13) with 

(pkiz) (pk+i{z) 
Mt) <Pk+i{t) ' 



{z - t)qk{z) 



and choose t so that (j)k{t) / for = 0, 1, . . . , A'o — 1. Instead of (2.14) consider the product 



No-i r 

n 

k=0 '- 



<Pk+l{t) 



k{t) XkXk+i 



n-l 

n 

k=No 



i>k(S>) XkXk+1. 



and take the limit t ^ so that t avoids the set fi. This proves equation (2.10) under the condition 
Dn{f) 7^ 0, n > 0. We extend the result to a weaker condition Dn{f) ^ 0, n = Nq, Nq + 1, . . . , and 
thus complete the proof of (2.10), by using the fact that Dn{f) 7^ 0, n = 1, 2, . . . , for positive / on 
the unit circle, and by a simple continuity argument in aj and Pj (cf. [28]). The case of z~^f{z) is 
dealt with similarly. □ 

Wc will now establish a connection between a Hankcl determinant with symbol on a finite interval 
and a Toeplitz determinant. First we need a theorem due to Szego on a relation between polynomials 
orthogonal on an interval of the real axis and those orthogonal on the unit circle. Szego considered 
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positive weights on the unit circle, but his theorem is transferred to the general case without much 
change: 

Lemma 2.5. Let f{z) have the property /(e*^) = /(e~'^), < < 27r and let 

w{xj = --7, x = cosO. 

I smu\ 

Assume that Dn{f) ^ 0, n > Nq, Nq > 0. Then the polynomials Pn{x) = ^Cnx"^ + • • • ; n = 
Nq^Nq + 1, . . . exist which are orthonormal w.r.t. weight w{x) on [—1,1], i.e., 

Pn{x)x"^w{x)dx = ^^^Snm, m = 0, 1, . . . , n, n > A^o, 

-1 

and, for n = Nq, Nq + 1, . . . , there hold the following expressions in terms of the polynomials (j)n{z) 
orthogonal w.r.t. f{z) on the unit circle: 



/: 



(2.15) x„ = 2"x2n 



1 — a2n-l 



27r 



(2.16) Pnix) = (2^)n(i_^^^_^) (^2n(-^^) + $L(-^)), n > Nq, 

where Pn{z) = Pn{z) / Xn, ^n{z) = (l)„{z)/Xn, Ki^) = ^"^n(^~^), ^n-l = -*n(0), n > Nq. 

Proof. The condition on Toeplitz determinants immediately implies the existence of the polynomials 

4'n{z) = Xn-z" + ' ■ ■ ; Xn 7^ 0) = -^0) -^0 + 1) • • • orthogonal w.r.t. f{z) on the unit circle. By Remark 
(b) above, in the present case of /(e*^) being an even function of 0, we have (pniz^^) = (pn{z~^) for 
all n > Nq. Now the proof is the same as the argument in the proof of Theorem 11.5 in [38], and 
we obtain 

(2.17) Pn{x) = -==^==(z-"</.2„(z) + z''<P2n{z-')). 

^y2^T[l - a2n-l) 

Note that 1 — a2n-i 0, n = Nq, Nq + 1, . . . as follows from (2.6) which in our case can be rewritten 
in the form 

Xn = (1 - <^n+l(0)/Xn+l) Xn+1- 

We now easily obtain the statement of the lemma from (2.17). □ 

Note that the recurrence relation (2.5) can be easily rewritten in terms of the monic polynomials 
in the form (for /(e*^) an even function of 9 and with z is replaced by z~^): 

(2.18) ^n+i{z) = z^n{z) - a„$;(z). 
Replacing here again z by z~^ and multiplying both sides by z""*"^ we obtain 

(2.19) K+li^) = Ki^)-(^nZ^n{z). 

Now we are ready to formulate and prove 

Theorem 2.6. [Connection between Toeplitz and Hankel determinants] Let Nq > and Dn{f) ^ 
for all n > Nq. Let the weights f{z) and w{x) be related as in Lemma 2.5. Let, moreover, 

Dn(w{x)) = det( f x^^'^w{x)dx\ , n = A^o, -^o + 1, . . . 
be the Hankel determinant with symbol w{x) on [—1, 1]. Then, with ^n{z) = 4'n{z)/xn, we have 

(2.20) D^[w[x)f = ^(1 + ^2n{0)?^^^§^y n = Nq, Nq + 1, . . . 
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Proof. Assume first A^o = 0. Take equation (2.16) with n = A; + l and apply the recurrence relations 
(2.18,2.19) with n = 2k + 1 to ^2k+2{z) and ^2fe+2(-^)' respectively. We then obtain 

Pfe+i(z) = (22)-^-i(z$2ik+i(^) + 

Now apply again the relations (2.18,2.19) with n = 2A; to ^2k+i{z) and ^2k+ii^) l^ere, respectively. 
The result can be written in the form 

^UZ) = ^^Pk-,1{X) - z(^^2kiz), 

1 - za2k 1 — za2k 

where we assume that z ^ and 1 — za2k / 0. On the other hand, from (2.16) with n = k 

$*fc(z) = (2z)'=(l - a2k-i)Pk{x) - ^2k{z) 

Equating the r.h.s. of the last two equations, wc obtain 

(2.21) (^2 - l)<^2k{z) = {2zf+^Pk+i{x) - {2zf{l - a2k-i){l - za2k)Pk{x). 

Setting here z = 1 (recall from the proof of Lemma 2.5 that 1 it a„ 7^ 0, n = 0, 1, . . . ), we obtain 

il-a2k-i){l-a2k) = 2^^. 

Note that it is a general property of orthogonal polynomials that Pk{z) and Pk+i{z) cannot have 
a zero in common. Similarly, setting z = —1, we have 

{l-a2k-i){l + a2k) = -2^^^^. 

The product of these two equations yields 

(1 - a2k-i) (1 - a2k) - -4- 



Pk{l) Pki-l 
By the relation (2.6), we can substitute here 

2 

l-a2ifc 



.2 ^ f X2k \ 



\X2k+lJ 



which gives 



X2k+l\ -Pfc+l(l)-Pfc+l(-l) 



(2-22) (l-a2fe-i)^ = -4, prnP^ n 

X2k J Pk[)-)Pk{-^) 

This equation together with (2.15) and the well-known expression for Dn{w{x)) in terms of the 
leading coefficients Xj^^ implies 

(2.23) DM = n-*-' = 4...-.|^(i)^(-i) n xf = ^jr^^^i^M 

Now using (2.16), we obtain 

(-1)" 
4"-i(l-a2„-i) 

and thus finish the proof for A'^o = 0. (Note that (2.22) implies that $2n(=tl) 7^ 0-) The extension 
to an arbitrary Nq > is carried out as in the proof of Lemma 2.4. □ 



P„(1)P„(-1) = 'J^ ^$2n(l)$2n(-l), 



We will also need a connection between Hankel and Toeplitz+Hankel determinants. We borrow 
the idea of the next statement from [39, 29, 2]. 
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Lemma 2.7. [Connection between Hankel and Toeplitz+Hankel determinants] Let fj be the Fourier 
coefficient fj = ^ /^f" f{e"^)e-'^'^de. Let f{e"^) = /(e"^^). Then, for n = 1,2, ... , 

(2.24) det(/,_fc + fj+k)]-kLo = DMx)), 

where Dn{v{x)) is the Hankel determinant with symbol v{x) = /(e*^*-^^)/\/l — x^, x = cos9 on 
[—1, 1]. Furthermore, again in terms of Hankel determinants with symbols on x & [—1, 1], 

(2.25) det(/,_fe - fj+k+2)]rkLo = —Dnifie''^^^)Vl^), 

(2.26) det(/,_, + fj+k+i)]-,Lo = ^^"(/(^''^'Vr^^' 



nn —n / 1 ™ 

(2.27) det(/,_, - /,+.+i);;,io = -;^^Me^'^^^)^^). 

Proof Since f{e'^) = /(e"*^), note that for j,k = 0,1,..., 

-| p27r -1 /•27r 

(2.28) - / f{e'')cosjecoskede = — / /(e^«)(e-^(^+'=) + e-^(^-^))d^ = /,_fe + /,+fc. 
^ Jo Jo 

Therefore, using the standard expansion (where only the first coefficient is needed to be known 
exphcitly) in non-negative powers of the cosine, 

(2.29) cos ke = 2^-^ cos'' 9 + Ck-2 cos^-^ 9 + Cfc_4 cos''-^ 9 + -- - , 
we obtain 

r2n \ n-1 



(2.30) det(/,_fc + fj+k)]Jilo = det (^^ ^ " f{e'^) cos j9 cos k9d9^ 

= 2i+2+- +"-2 det (^^ j'^J f{e'^) cos i^cos'^ 9d9^ 
= 2("-i)("-2) det (^^ j^J f{e'^) cos^' ^cos'' 



j,k=0 

9d9 



j,k=0 

n / \ n—1 



= 2("-i)("-2) (J-^ det f{e'^) cos^ 9 cos*^ 9de 



j,k=0 



Changing the variable x = cos 9, d9 = -dx/yjl — x^, we immediately obtain 

r)-n?-2n+2 / i-l \ f(JS(x)\ 



2n^-2n+2 /pi \ 

(2.31) det(/,_fe + fj+k)]-kLo = det i^J^^ v{x)x^+''dxj , vix) 



which is (2.24). 

Similarly, using the observations 

r2TT 



(2.32) - / " f{e'^) sm{j + 1)9 sm{k + l)9d9 = fj^k - fj+k+2, 

Jo 

(2.33) - r f{e'') cos(j + 1/2)^ cos(A; + l/2)9d9 = fj_k + fj+k+u 

Jo 

(2.34) - / f{e'^) sm{j + 1/2)9 sin(A; + l/2)9d9 = fj^k - fj+k+i, 

TT Jo 
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and the expansions in non-negative powers of the cosine of the quantities 

sin(A; + l)6' cos(A; + 1/2)6' sin(A; + 1/2)6* 
sin 6* ' cosf ' sin I 

we obtain (2.25), (2.26), and (2.27). □ 

Finally, we list some properties of Barnes' G-function (see [3, 40]) we need below. The G-function 
is an entire function defined, e.g., by the product: 

(2.35) G(z + 1) = (27r)^/2e-(^+i)^/2-^^^V2 -Q (^1 + |j ^-zWl{2k)^ ^ C 

fc=i 

where is Euler's constant. G(z) satisfies the recurrence relation: 

(2.36) G(z + 1) = r(z)G(z), G(l) = l, 
where r(z) is Euler's G-function. The following representation is useful 

(2.37) / lnr(x + l)dar = Jln27r-^^^-til + zlnr(z + l)-lnG(z + l). 
Jo 2 2 

There holds the identity: 

(2.38) 21nG(l/2) = (1/12) In 2 - In Vtt + 3C'(-1), 

where C^{x) is the derivative of Riemann's ^-function. We will also need a doubling formula given 

by 

(2.39) G(2z)7r^G(l/2)2 = G{zfG{z + l/2)2r(z)2(2^-i)(^-i) . 

3. RiEMANN-HlLBERT PROBLEM 

In this section we formulate a Riemann-Hilbert problem (RHP) for the polynomials ^^(-z), (pki^)- 
We use this RHP in section 5 to find asymptotics of the polynomials. 

Let the weight f{z) be given on the unit circle (which, oriented in the positive direction, we 
denote C) by (1.2). Suppose that the system of orthonormal polynomials satisfying (1.16) exists. 
Consider the following 2x2 matrix valued function Y^'^^z) = Y{z): 



(3.1) yW(z) = 



It is easy to verify that Y{z) solves the following Riemann-Hilbert problem: 

(a) Y{z) is analytic for z € C \ G. 

(b) Let z e G\U"io'2i- ^ has continuous boundary values Yj^[z) as z approaches the unit circle 
from the inside, and Y-{z), from the outside, related by the jump condition 

(3.2) Y+{z)=Y^{z)(\ ^"Y^^A, zGG\U-o^,-. 



(c) Y{z) has the following asymptotic behavior at infinity: 

z^ ( 
Z J J \^ z~ 



(3.3) Y{z)=[l + 0[-\](i ° ), asz^oo. 
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(d) As z ^ Zj, j = 0,1, ... ,171, z e C\C, 



(3.4) 



Yiz) = (; 




:) 



if aj 7^ 0, 



and 



(3.5) 



^/0(1) 0{\n\z-z,\)\ 
VO(l) 0{ln\z-zj\)J^ 




0, Pj ^ 0. 



(Here and below 0(a) stands for 0(|a|).) 

A general fact that orthogonal polynomials can be so represented as a solution of a Riemann- 
Hilbcrt problem was noticed in [24] (for polynomials on the line) and extended for polynomials on 
the circle in [1]. This is important because it turns out that the RHP can be efliciently analyzed for 
large by a steepest-descent-type method found in [17] and developed further in many subsequent 
works. Thus, we first find the solution to the problem (a)-(d) for large k (applying this method) 
and then interpret it as the asymptotics of the orthogonal polynomials by (3.1). 

The solution to the RHP (a)-(d) is unique. Note first that dety(z) = 1. Indeed, from the 
conditions on Y{z), dety(z) is analytic across the unit circle, has all singularities removable, and 
tends to 1 as z ^ oo. It is then identically 1 by Liouville's theorem. Now if there is another 
solution Y{^z\ we easily obtain by Liouville's theorem that Y{z)Y{z)~^ = 1. 



In this section we construct an asymptotic solution to the Ricmann-Hilbert problem (a) (d) of 
Section 3 for large A; = n by the steepest descent method. All the steps of the analysis are standard 
apart from construction of the local parametrix near the points Zj. We always assume that f{z) is 
given by (1.2). In this section we also assume for simplicity that zq = 1 is a singularity. However, 
the results trivially extend to the case uq = Pq = 0. In this section and the next one, we further 
assume that V{z) is analytic in a neighborhood of the unit circle. 

The first step is the following transformation, which normalizes the problem at infinity: 



4. Asymptotic analysis of the Riemann-Hilbert problem 



(4.1) 




Prom the RHP for Y{z), we obtain the following problem for T{z): 

(a) T{z) is analytic for z G C \ C 

(b) The boundary values of T(z) are related by the jump condition 




(c) T{z) = 1 + 0{l/z) as z^oo, 
and the condition (d) remains unchanged. 

Now split the contour as shown in Figure 1. Define a new transformation as follows: 

'T{z), for z outside the lenses, 



(4.3) 





for |z| > 1 and inside the lenses. 
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Figure 1. Contour for the S'-Riemann-Hilbert problem (m = 2). 



(b) The boundary values of S{z) are related by the jump condition 
where the minus sign in the exponent is on T,j, and plus, on T,", 

(c) S{z) = 1 + 0{l/z) as z^oo, 

(d) As z ^ Zj, j = 0, ... ,m, z E C\C outside the lenses, 



(4.4) Siz) = 
if aj ^ 0, and 

(4.5) S{z) 



0(1) 0(1) +0(|Z- Zjf"^' 

0(1) 0{\)+0{\z- 



_ fO{l) 0{ln\ 
~ \Oil) 0(ln| 



Z — Zi 



if aj = 0, Pj 7^ 0. The behavior of S{z) for z —>■ zj in other sectors is obtained from these 
expressions by application of the appropriate jump conditions. 

Let us encircle each of the points zj by a sufficiently small disc, 
(4.6) U,. = {z:\z-Zj\<s}, 

We see that, outside the neighborhoods Uzj, the jump matrix on Sj, j = 0, . . . , m is uniformly 

exponentially close to the identity. We will now construct the parametrices in C \ (U^g^Zj) 
■ We match them on the boundaries dUz- , which yields the desired asymptotics. 

4.1. Parametrix outside the points Zj. We expect the following problem for the parametrix N 

inC\U™o^.r 

(a) N[z) is analytic for z G C \ C, 
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(b) with the jump on C 



N.{z) 





-/(^) 







^ G C\U 



(c) and the following behavior at infinity 



N{z) 



i + o[- 

z 



as z 



GO. 



One can easily check directly that the solution to this RHP is given by the formula 

'V{zY\ \z\ > 1 

' l\ 



(4.7) 

where the Szego function 
(4.8) 



Niz) 



V{z) 



-1 0, 



T>[z) = exp 



-/ 

2m Jc 



kl < 1 ' 



ds, 



Ic s-z 

is analytic outside the unit circle with boundary values satisfying ^^{z) = V^{z)f{z), z e C \ 

In what follows, we will need a more explicit formula for 'D{z). Calculation of the integral (with 
the help of (4.13) below) gives: 
(4.9) 

V{z) 

and 
(4.10) 

Viz) = 



exp 



2mJcs-z V -Zfce*'' / +^ ' LL 



\z\ < 1. 



exp 



^ / ^^^1 n 



Zk 



-CCk+Pk 



=6_(z)-in 

fc=0 



Zk 



-O!k+0k 



\z\ > 1, 



where Vq, h±{z) are defined in (1.8). Note that the branch of (z — 

^^)±afe+/3fc in (4.9,4.10) is taken 

as discussed after equation (4.13) below. In (4.10) for any k, the cut of the root ^^"fc+^fc is the line 
9 = 9k from z = to infinity, and 9k < arg z < 2tt + 9k. 

4.2. Parametrix at Zj. Let us now construct the parametrix Pzj{z) in Uzy The construction is 
the same for all j = 0, 1, . . . . We look for an analytic matrix-valued function in a neighborhood of 
which satisfies the same jump conditions as S{z) on S H f/^^. , the same conditions (4.4,4.5) as 
z — > Zj, and, instead of a condition at infinity, satisfies the matching condition 

(4.11) Pz,{z)N-\z) = I + oil) 

uniformly on the boundary dUz^ as n — ^ oo. 
First, set 



(4.12) 



C = nln— , 



where Inx > for x > 1, and has a cut on the negative half of the real axis. Under this transforma- 
tion the neighborhood Uzj is mapped into a neighborhood of zero in the ^-plane. Note that C,iz) is 
analytic, one-to-one, and it takes an arc of the unit circle to an interval of the imaginary axis. Let 
us now choose the exact form of the cuts S in Uz^ so that their images under the mapping C,iz) are 
straight lines (Figure 2). We add one more jump contour to S in Uz^ which is the pre-image of the 
real line Fa and Fy in the ^-plane. This will be needed below because of the non-analyticity of the 
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function \z~ Zj\ . Note that we can construct two different analytic continuations of this function 
off the unit circle to the pre-images of the upper and lower half ^-plane, respectively. Namely, write 
for z on the unit circle, 

(4.13) h^,{z) = \z- zjl-^ = {z- zj^/'iz-' - zj'r/' = z = 

where £j is found from the condition that the argument of the above function is zero on the unit 
circle. Let us fix the cut of {z — zj)"^ going along the line 9 = 6j from zj to infinity. Fix the branch 
by the condition that on the line going from zj to the right parallel to the real axis, 8iTg(z — Zj) = 27r. 
For in the denominator, < argz < 27r (the same convention for roots of z is adopted in 

(4.15,4.17) below). Then, a simple consideration of triangles shows that 



(4.14) Ij = 



37r, < 6' < 9j 



Thus (4.13) is continued analytically to neighborhoods of the arcs < 9 < 9j, and 9j < 9 < 27r. 
In Uzj, we extend these neighborhoods to the pre-images of the lower and upper half C-plane 
(intersected with ({Uzj)), respectively. The cut of haj is along the contours Fa and Fy in the 
C-plane. 

For z — Zj, C = n{z — Zj)/zj + 0{{z — Zj)^). We have < arg^ < 27r, which follows from the 
choice of arg(z — Zj) in (4.13). 

We now introduce the following auxiliary function. First, for j ^ 0, 



V(z) f 7 \ 

(4.15) F,{z) = e^J[[-] \{K,{z)9p,{z)^'^ 
k=o ^^''^ k^j 



i -inaj C e I II V VI 
gi^a, , ( e III, IV, VII, VIII ' ^ 

The functions gp^{z) are defined in (1.4). The case of Uz^ is slightly different because of the branch 
cut of z^'^ and z"'^ going along the positive real half-line. Let a step function 

^/i ifiA - f \ Je-^-*, argz>0 



and define 



e*^*, argz<27r' 

V(z) / r \ ^fe/^ 

(4.17) Fo(z) = e^n r \{h^.i^)mM'^ 

u-c \Zk/ , ,„ 



„-i7rao^ (el, 1 1 

e^-(°o-/3o)^ (^e III, IV 

c^vii,viii 

It is easy to verify that Fj{z), j = 0,1, . . . is analytic in the intersection of each quarter ^-plane 
with C{Uzj) and has the following jumps: 

(4.18) F,,+ (z) = F,,_(z)e-2™.- C G Ti; 

(4.19) F,,+ (z) = F,,_(z)e2--^ C G Tg; 

(4.20) F,-+(z) = F,-_(z)e""^ C G U Fr. 
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Figure 2. The auxiliary contour for the parametrix at Zj 



Comparing (1.2) and (4.15), and using the analytic continuation (see (4.13)) for /(z) off the arcs 
between the singularities, we obtain the following relations between /(z) and F^(z): 

(4.21) Fj{zf = /(z)e"2™--i' 

(4.22) 



F,{zf = /(z)e2-"^5^;(z) 



C G III, IV, VII, VIII. 



for j 7^ 0. If j = the same relations hold with the functions g^^^ {z) replaced by g^^ (z) . 
We look for P^- (z) in the form 

(4.23) P,. {z) = E{z)P^^^ (z)Fj (z)-^3^±"^3/2^ 

where plus sign is taken for \z\ < 1 (this corresponds to ^ G I, II, III, IV), and minus, for \z\ > 1 
(C G V, VI, VII, VIII). The matrix E{z) is analytic and invertible in the neighborhood of Uzj, 
and therefore does not affect the jump and analyticity conditions. It is chosen so that the matching 
condition is satisfied. 

It is easy to verify (recall that Pzj{z) has the same jumps as S{z)) that P^^\z) satisfies jump 
conditions with constant jump matrices. Set 

(4.24) pW(z) = *,(C). 
Then ^j{Q) satisfies a RHP on the contour given in Figure 2: 



(a) is analytic for C G C \ ^]=iTj. 

(b) satisfies the following jump conditions: 



(4.25) 

(4.26) 
(4.27) 
(4.28) 



*.,+(C) 





,■+(0 = *,,-(C)(_A/3. 

0-3 



-j7r/3j 



*,-+(C) = *,-,-(C)e^""^- 





gj7r/3 



for C G Fa U Fr 



for C G Fi, 
for C G Fs, 
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for C € r2 with plus sign in the exponent, for G with minus sign, 

(4.29) *,-,+(C) = *,-,-(C) (^^±inil+2a,) 1 

for C £ Ts with plus sign in the exponent, for ^ G Fg, with minus sign, 
(c) As C ^ 0, C e C \ U^=irj outside the lenses, 

/0(O) 0(O) + 0(C"°^ 



(4.30) ^j{z) 

if 7^ 0, and 

(4 31) *(.)-f^^^^ ^^^"1^1)^ 

^^■-^^^ ^^^^^^-1^0(1) 0(ln|C|); 

if aj = 0, Pj ^ 0. The behavior of '^j{z) for — in other sectors is obtained from these 
expressions by application of the appropriate jump conditions. 

Wc will solve this problem explicitly in terms of the confluent hypergcomctric function, i/j{a,, c: z) 
with the parameters a, c determined by a^, [3j. A standard theory of the confluent hyper geometric 
function is presented, e.g., in the appendix of [28]. 

Denote by Roman numerals the sectors between the cuts in Figure 2. The following statement 
holds. 

Proposition 4.1. Let aj ± f3j / —1, —2, . . . for all j. Then a solution to the above RHP (a)-(c) 
for ^j(C); < argC < 27r, is given by the following function in the sector I: 



(4.32) ^jiO = (0 = _ + 1 _ 2a„ C)e'-(/^^--3".0e-C/2 ra+j^y.) 

where ip{a,b,x) is the confluent hypergeometric function, and F(x) is Euler's T-function. The 
solution in the other sectors is given by successive application of the jump conditions (4-25-4-29) 
to (4.32). 

Remark 4.2. The functions C^"^ ) V'('j^) Oj s-i^d ip{a, b, e~^'^Q are defined on the universal covering 
of the punctured plane C, € C\ {0}. Recall that the branches are fixed by the condition < arg^ < 
27r. 

Proof. The condition (c) is verified in the sector / by applying to (4.32) the standard expansion of 
the confluent hypergeometric function at zero (see, e.g., [11]), namely, 

(4.33) V(o, c, x) = (1 + 0{x)) + (1 + ' 

X — ^ 0, c ^ Z, 

or, to cover also the integer values of c: 

'^^^x^-''{l + 0{x\nx)) + 0{l), 3f?c>l 
^^i^ (1 + Oix)) + ^^x^-'^ (1 + Oix)) , 3fic = 1, c 7^ 1 



(4.34) tp{a, c, x) 



r(l+a-c) 

-^(lnx + ^-2-fE)+0{x\nx), c=l 
^ ra+I-c) (1 + 0{x\nx) + 0(x^-^)) , 3fic < 1 



0, 
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where = 0.5772 ... is Euler's constant. 

We verify the condition (c) similarly in the other sectors. 

To verify (6), reduce the contour of Figure 2 to the real line, oriented from right to left, by 
extending the sectors I and IV and collapsing the jump conditions. We then obtain the following 
reduced RHP: 

(4.35) <;'(C) = •i':^(QWi' = K) {^^^^_ „.), , < < 0; 

*Sf]:'(c) = *<:!(OA-'Js-'^f' JaJ. = *<:!(c) ('-";*'-°'' ^''^■Jt.?'^ . < > «. 

where the jump matrices Jk correspond to jumps on the contours F^, k = 1, ... ,8 as defined in 
(4.25-4.29). 

The confluent hyper geometric function possesses the following transformation property on the 
universal covering of the punctured plane: 

(4.36) c, e-2-C) = e^-^a, c, C) - r(a)F(a - c + 1) "'"""'^'^^^ ~ 

This property is proved in the appendix of [28] (equation (7.30)). 

Taking ^'f^(C) given by (4.32) and applying to it the jump condition for C < 0, we obtain using 

(4.36) and the standard properties of F-function the following expressions for the first column of 

(4.37) ^^;j^{Q = OVK+/3,-,l + 2a„e-2-C)e-^/2 

(4.38) *(.2)(C) = - aj + Pj, 1 - 2aj, e-^-^C)e-^-^-e-</^ ^^^ + ^ 

i [aj pj) 



The second column is 



(4.39) (C) = ^fj,i0e-'-"^, ^fj,\0 = *S2(C)e-^'^"^- 

Now applying to this function the jump condition for C > and using again (4.36), we obtain 
(note that as a result of these manipulations we moved C — e^'^X) 

(4.40) *f)(C) = *,(C) 

with < argC < t^, i-e. the ^^'^''(C) we started with. Thus, (4.32, 4.37) is a solution to the 
reduced RHP given by the jump condition (4.35). Therefore, (4.32, 4.37-4.39) give a solution to 
the original RHP for * in the sectors / and IV, respectively; and the solution in the other sectors 
is reconstructed using (4.25-4.29). Proposition 4.1 is proved. □ 

We will now match this solution with N{z) on the boundary dUzj for large n. The limit n — oo, 

z G dUzj, corresponds to ^ ^ oo, therefore we need the asymptotic expansion of ^j(C)- We use 
the classical result (e.g., [11] or Eq.(7.2) of [28]) for the confluent hypergeometric function: 

(4.41) '0(a,c,x) =x~"[l-a(l + a-c)a;~^ + 0(x~^)], ^ oo, -37r/2 < argx < 37r/2. 

Note that these asymptotics can be taken both for il){a, c, C,) and '4^{a, c, e~'^^C,) for C, ^ I. We apply 
this result to (4.32) and thus obtain the asymptotics of the solution in the sector I. The "proper" 
triangular structure of the jump matrices implies that these asymptotics remain the same in the 
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sector II as well, namely: 

Furthermore, applying the jump matrices, we obtain the following asymptotics for ^'j(C) in the other 
sectors (here ^'^■'(C) stands for the analytic continuation of the r.h.s. of (4.42) to < arg(^ < 27r) 
as C — >■ oo: 

(4.43) ^f'\Q = ^f\Q = vI//)(C)e^™^'^«, 

(4.44) ^rW(0 = *f^)(C) = *f(C)(^-!l;3, -^^'^^^^e--'^^-^^ 

(4.45) *r)(c)=*r"^(c) = *f(c)(A- "'o'^'O- 

Now substituting these asymptotics into the condition on E: 

(4.46) P^.{z)N-^{z) = E{z)^j{()Fj{z)-''^z^'"'^/^N-\z) = I + o(l), 
we obtain 

(4.47) E(z) = iV(.)C^.-3F- (z)z7'^'^^/^ Q ' for C G I, /I, 

(4.48) Eiz) = Ar(^)c/^.-3F-3(,)^-n.3/2 e ^ ^,.(^'+3a,)J - for C e m, /y, 

(4.49) £;(^) = iV(.)C-^^-3j.^^3(,)^-3/2|'^^_^J^^^^^^^ e ^ forCGF,V^I, 

(4.50) i?(z) = iV(z)C-^^-3i.-3 (^)^-3/2 ^ y for C G VII, VIII. 

The dependence on z enters into these expressions only via the combination V{z) / {(^^^ Fj(^z)) 
for \z\ < 1 (i.e., ^ G I, II, III, IV) and the combination D (z)Fj(z)/C'^J- for |z| > 1 (i.e., C e 
y, 1//, y//, VIII). Expanding the logarithm in (4.12) in powers of u = z — Zj, we see immediately 
from (4.9,4.10,4.15,4.17) that the mentioned combinations, and therefore E{z) have no singularity 
at Zj. Thus E{z) is an analytic function in Uzy In what follows, we will need more detailed 
information about the behaviour of some of these combinations as — 0. Namely, it is easy to 
obtain from (4.12,4.9,4.15,4.17) and (4.13) that 

(4.51) Fj(z) = ?7je-3*™^/2z7"^^x"^(l + 0(u)), u = z-Zj, (el, 



where 



(4.52) r/j- = e^(^^)/^exp< 



.k=0 k=o+l I I k^j 
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and 

(4.53) ix^^y = /^^'"^"^-'^^^"-'^^(l + 0{n)), u = z-zj, (el, 

- = i^-m) {-^ (I- - i.-) ) n (1)""^ - ^'1^'- 

To derive (4.54), we used, in particular, the factorization (1.8). 

It is seen directly from (4.47-4.50) that det£;(z) = e'''^°'^'^^\ Note that as follows by Liouville's 
theorem from the RHP, det ^'j(C) = e~^'^^'^^~^^^: this function has no jumps, the singularity at zero 
is removable as ^aj > —1/2, and the constant value follows from the asymptotics (4.42). Combining 
these results, we see from (4.23) that det Pzj{z) = 1. Comparing the conditions (4.30,4.31) and 
(4.4,4.5), we see that the singularity of S{z)Pz- (z)^^ at z = zj is at most 0{\z — Zjp"^) or 0(ln \z — 
Zjl'^). However, by construction of P^^, the function S{z)Pzj{z)~'^ has no jumps in a neighborhood 
of Uz- and hence this singularity is removable. Thus, S{z)Pz {z)~^ is analytic in a neighborhood 
oft/.,. 

Note that the error term in (4.46) o(l) = n"^'3j'^;^0(n"i)n^'3^'^^ It is o(l) for -1/2 < < 1/2. 

This completes the construction of the parametrix at zj: it is given by the formulae (4.23,4.24,4.47- 
4.50) and Proposition 4.1. 

Considering further terms in (4.42), we can extend (4.46) into the full asymptotic series in inverse 
powers of n. For our calculations we need to know explicitly the first correction term: 



(4.55) Pz^{z)N-^{z) =1 + Ai(z) + n-^^^'^30(i/j^2)^SR/3,a3^ 



Ai(z) 



-(af-/3„ 



r(l+aj+/3,) 



V r(a,+/3,) W^3Fiiz) I 



zedz{I), aj±Pj^ -1,-2,. 



where dz{I) is the part of dUz^ whose C-image is in /. As a consideration of the other sectors 
shows, this expression for Ai(z) extends by analytic continuation to the whole boundary dUz-- As 
follows from (4.53), it gives a meromorphic function in a neighborhood of t/. with a simple pole 



clt Z — Zj . 



The error term 0(l/n ) in (4.55) is uniform in z on dUz ■ 



4.3. R-RHP. Throughout this section we assume that aj±f3j / —1, —2, ... for all j = 0, 1, . . . , m. 

Let 

(4.56) R(z) = l^^'^^''^'^' ^eC/oo\S, C/oo = C\U-oC/.,, 

zeUzjX'S, j = 0,...,m. 

It is easy to verify that this function has jumps only on dUzj , and parts of , lying outside 
the neighborhoods Uzj (we denote these parts without the end-points S""*). The contour is shown 
in Figure 3. Outside of it, as a standard argument shows, R{z) is analytic. Moreover, we have: 
R{z) = I + 0{l/z) as z^ oo. 
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The jumps of R{z) are as follows: 

(4.57) R+{z) = i?_(z)iV(z) (^^^^_\^„„ flN{zr\ 2;GS-\ 

(4.58) RAz) = i?-(^)^(^) (;(^)-i^n j)iV(z)-\ 

(4.59) R+{z) = R_{z)Pz.{z)N{z)~^, z G St/^. \ {intersection points}, 
J = 0, . . . , m. 

The jump matrix on E""^*, S can be estimated uniformly in aj, (5j as / + 0(exp(— en)), where 
£ is a positive constant. The jump matrices on dU^^ admit a uniform expansion in the inverse 
powers of n conjugated by n'^j'^3^~"'^3/2 ^^-^^ g^^^ term is given explicitly by (4.55)): 

(4.60) 7 + Ai(z) + A2(^) + --- + Afe(z) + A2„ z&dU,.. 
Every Ap(z), aJ,''^(2;), p = 1, 2, . . . , z G Uj'^odU^j is of the form 

m 

(4.61) ^a7'^^0(n-naf , a,- = n'^^z/"/', 

j=0 

it is of order n'^^^^j l^l^jhp. 

To obtain a standard solution of the /2-RHP in terms of a Neumann series (see, e.g., [18]) we 
must have n'^^'^^J I^Rfthi = o{l), that is 3f?/3j G (-1/2, 1/2) for all j = 0, 1, . . . ,m. However, it is 
possible to obtain the solution in the whole half-closed interval 3fi/3j G (—1/2, 1/2], j = 0, 1, . . . , m, 
and moreover, in any half-closed interval of length 1. 

Consider the transformation 

(4.62) R{z) = n'^'''^R{z)n-^''^, 
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where 



^ (miiij- + (9 + 5)) > if several /Sj ^ 0, and ^/Sj G (g - 1/2, g + 1/2], g G 



(4.63) LV = < 



i [maxj + - 5)) , if several /3j 7^ 0, and ^/3j G [g - 1/2, g + 1/2), g G 

if there is only one nonzero Pj^ 
0, ifan/3j=0. 



will "shift" all 3f?/?j inside the interval (—1/2,1/2). Recall that = 1 is not considered if both 

ao = and /3o = 0. 

Now in the RHP for R{z), the condition at infinity and the uniform exponential estimate / + 
0(exp(— en)) (with different e) of the jump matrices on S is preserved, while the jump 

matrices on dUz^ have the form: 

(4.64) / + n'^''^Ai{z)n-'"'' + ■■■ + n'^"^ Ak{z)n-'"'^ + n'"'' A^^l^{z)n-'"'' , z G dU^., 
where the order of each rf"'' A^{z)n-'^'''\ rt^"'' A^y\z)n-'^''^ , p = 1, 2, . . . , z G WJLodU^. is 

This implies that the standard analysis can be applied to the i?-RHP problem in the range G 
{q-l/2,q + 1/2], j = 0, 1, . . . , m (or g - 1/2 < ^pj < q + 1/2, j = 0,1, . . . m) for any g G R, and 
we obtain the asymptotic expansion 

k 

(4.65) Riz)=I + J2Rpi^)+Rtlii^)^ k = l,2,... 

p=i 

In our case the error term 

(4.66) R^j^liz) = 0{\Rk+i{z)\) + Oi\Rk+2{z)\). 

The functions Rj{z) are computed recursively. In this paper, we will need explicit expressions 
only for the first two. Accordingly, set k = 2. The function Ri{z) is found from the conditions that 
it is analytic outside dU = U^^dUzj, Ri{z) ^ as z ^ 00, and 

(4.67) Ri,+ {z) = Ri,-{z) + n^''^Ai{z)n-'^''^, z G dU. 
The solution is easily written. First denote 

(4.68) i?p(z) = n-'^^3^p(z)n'"'^^ R^;\z) = n''"'^ R^;\z)n'"'\ 
and write for R: 

Ai{x)dx 



(4.69) R,i^) = 4-L 



du x-z 

"lELo54t^-Ai(.), zeUz,, i = 0,l,...,m.' ^^-^.-0^^^- 

where the contours in the integral are traversed in the negative direction, and A/, are the coefficients 
in the Laurent expansion of Ai (z) : 

(4.70) Ai{z) = ^^ + Bk + 0{z-Zk), z^zk, k = 0,l,...,m. 

Z Z}^ 
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The coefficients are easy to write using (4.55) and (4.53): 



(4 71) A.-A<»>-i( 

An expression for Bk is also easy to find, but it is not needed below. 

The function i?2 is now found from the conditions that i?2(-2) — as 2; — oo, is analytic outside 

dU , and 

(4.72) ^2,+ (-z) = ^2,-(^) +-Ri-(2)n'^'''Ai(z)ra-'^''3 +„'^^3^2(^)ra"'^''^ z G dU. 
The solution to this RHP is 

(4.73) R2{z) = ^ f (Ri -{x)n''^'Ai{x)n-'"'' + n""^' A2{x)n-''^A 

27rz Jqu \ ' / X- z 

Further standard analysis (cf. (4.66)) shows that the error term 

0{5/n) + 0(^2) o(s maxk 



(4.74) 4\z) 



O (<5 maxfc 0{d/n) + 0{d^) 



where S is given by (1.20). 

In particular, as is clear from the above, if there is only one nonzero Pj^, we obtain the expansion 
of R{z) purely in inverse integer powers of n valid in fact for all G C, aj^ ± Pj^ ^ —2, . . . 

It is clear from the construction and the properties of the asymptotic series of the confluent 
hypergeometric function that the error terms i?^!^''(z) (4.66), and in particular (4.74), are uniform 
for all z and for fij in bounded sets of the strip (7 — 1/2 < < (7 + 1/2, j = 0, 1, ... m, (or 
q — 1/2 < 3ft/3j < g + 1/2, j = 0, 1, . . . m) for ctj in bounded sets of the half-plane Jfta^ > —1/2, 
and for ctj it /3j away from neighborhoods of the negative integers. Moreover, the series (4.65) is 
differentiable in a^, 

For future use note that if V{z) = Vr{z) + {V{z) — Vr{z))h, h G [0, 1], and Vr{z) is analytic in a 
neighborhood of the unit circle, then the error terms are uniform in the parameter h G [0, 1] . 

5. Orthogonal polynomials. Proof of Theorem 1.8 

Using results of the previous section, we can provide a complete asymptotic analysis of the 
polynomials orthogonal with weight (1.2) on the unit circle with analytic V{z). In this section we 
will find the asymptotic expressions for Xn, 0n(O), and 0n(O)- 

First, it follows immediately from (3.1) that 

(5.1) xl-i = ~Yi^(0). 

Tracing back the transformations R^S^T^Y,we obtain for z inside the unit circle and 
outside the lenses: 

(5.2) Y{z) = T{z) = S{z) = R{z)N{z) = n-^''--'R{z)n'^''m{z) 

= n-'"'^[I + Ri{z) + R2{z) + R^^\z)]n'"'^N{z) 

= 'l + R,{z) + R2{z) + Rt\z)\ V{zr J) . 

Taking the 21 matrix element and setting z = we obtain 

(5.3) Xn-l = -^"^(0) = T){0)-' [1 + i?i,22(0) + R2,22m + 0(6 /u + 6^)] , 
where we used the estimate (4.74) for r'^\z). 
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-IX 



-Vo 



By (4.9) 

(5.4) p(0)-^ = exp 

Using (4.69) and (4.71) we obtain 

(5-5) RiMo) = - E ^ = E ("i - (^l) ■ 



Zk 



k=o fe=o 

Conjugating (4.73) with n'^'^^, setting there z = 0, and applying (4.69), we obtain: 



(5.6) 



2l3j 



From (4.71), 

(5.7) {AkAj)22 = ZkZj [{a) - p]){al - (3l)n-^ 



_^2(A-/3,-i) (z^y r(l + aj + /3,-)r(l + ak- (3k) f^] 



where are defined in (4.54). 

Substituting the last 3 equations into (5.3), we finally obtain (1.21). 
We now turn our attention to (?!>n(0). Using (3.1), we have 

(5.8) </>„(0) = XnY^l\Q) = Xn (uiOmor o))^^ = -Xn^(O)"' (^1.12(0) + . 

By (4.69,4.71), 



(5.9) 



Ak. 



12 



k=0 



Zk 



-. m 



3=0 



„-2/3,.n r(l + ai+/3j) 2 



and, recalling (4.66), we obtain (1.23). 

Finally, starting again with (3.1), we have 



(5.10) <f>n-l{0) 

We have 
(5.11) 



1- ^r'(^) 



Xn-l Z 



n— 1 



^ lim -l^(i2(z)P(z)-3z«-3)2i 



Xn-l '^^"O Z 
1 



Xn-l 



lim 2;i?i,2i(-2) + O 



n 



n 



max 

k n 



lim zR2iiz) = V Afe,2i = — E 



A;=0 



j=0 



i r(i + a,-/3,) 2 



and therefore, recalling (1.21), obtain (1.24). 

Note that uniformity and differentiability properties of the asymptotic series of Theorem 1.8 
follow from those of the -R-expansion of the previous section. 
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6. TOEPLITZ DETERMINANTS. PROOF OF THEOREM 1.13 



6.1. The case of analytic V{z). First, let V{z) be analytic in a neighborhood of the unit circle. 

(r) 

Consider the set /3j ' constructed in the proof of Lemma 1.12. We have to consider only the second 
class, i.e. = 1- We then have, relabeling /3^^^ according to increasing real part. 



,(0 



(r) 



3(0 



(6.1) - ■■■ - jv,yp ■ \ ^>-^^p+i - • • • - -^^i^m'-l ^ ^"-/^m'-^+l 



for some p, ^ > 0. Here w! is the number of singularities: m' = m + 1 if z = 1 is a singularity, 
otherwise m' = m. Now consider the symbol (not a FH-representation of /) / of type (1.2) with 
beta-parameters denoted by /3 and given by /3j = fi^p for j = 1, . . . , m' — and /3j = fi^ — 1 for 
J = m' — £ + 1, . . . , m'. It is easy to see that the original symbol / has (^^^) FH-representations 

in M. obtained by shifting any i out of £ + p parameters fij , say , . . . , , with the smallest real 
part to the right by 1. Thus, 



(6.2) 



/(^)=(-i)'n^. 



'^•xz-i...^-V/(z), 



j=0 



for appropriate Lj. 

Let us now and until the end of this section relabel /3j, ctj, Lj, and Zj according to increasing 
real part of Thus, in particular, 

(6.3) 3f?^i = . . . = 3f?^,+j, < 3f?^,+p+i . 

Assume that the set of all the minimizing FH-representations M. is non-degenerate (see Introduc- 
tion). This implies that ± /3j 7^ — 1, — 2, 

We now apply Lemma 2.4 (equation (2.9)) to finish the proof of Theorem 1.13. We need to 
evaluate the determinant n > Nq for a sufficiently large Nq > 0. First, from (3.1), tracing 
back the transformations of the RH problem and using (4.69,4.71) we obtain (cf. (5.8,5.9)) for the 
polynomials orthonormal with weight f{z): 



A 



in) 
fc,12 



(6.4) Mz)/Xn = V{z)-'pn{z), pn{z) = -J^ + O 



k=0 



Zk 



n 



n 



-2SR/3i-l 



This expansion is uniform and differentiable in a neighborhood of zero. A simple algebra shows 
that in the determinant 

(f>n+l{0)/Xn+l 
^( 



(6.5) 



0n(O)/Xn 



£MO)/Xr. 



£<f>n+l{0)/Xn+l 



(f>n+£-l(S>)/Xn+e-l 



£(f)n+e-l{0)/Xn+e-l 



^I^0n(O)/Xn ^FrT0n+l(O)/Xn+l ••• ^I^^n+^-l (0)/Xn+^-l 

all the terms with the derivatives of ^{z) drop out, and we have 



(6.6) 



p(o)- 



^rrTPn(O) J^rn-Pn+l (0) ••• J^Pn+l-li^) 

It is a crucial fact that the size C. of this determinant is less than the number of terms, i+p, in the 
expansion of (t)n{^)/Xn of the same largest order 0{n~'^^f^^~^) (see (1.23) with (5j replaced by (5j). 
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As |3fi/3j — 3f?/3jk| < 1, and aj it /3j 7^ —1, —2, . . . , j,k = 1, . . . , m', we obtain for the p'th derivative 
of p{z) from (6.4), (4.71), and (4.54) with P replaced by P 



(6.7) —p^+,{0) = s\J2^^ + O 



k=0 



n 



n 



-20?/3i-l 



n 



n 



-2g?/3i-i 



where 
(6.8) 



r(a,--/3,) 



■ exp < 



-ZTT 



i-1 m \ j 

fc=0 fc=J+l / I fcT^J 



Oik 



|2/3fc 



Substituting these expressions into the determinant we obtain 

(6.9) F„ = P(0)-^n^' E dA---d,,zl---z^,;'^^ n + 

=p(o)-^n^' E (inrfi,---di,(zi,---zi,)- n 

s=0 l<ii<i2<---<i^<£+p l<j<fe<^ 

Therefore, by (2.9), 

(6.10) I)„(^V(^)) = ^^/i ," ^n(/(^)) 

= (-1)"^P(0)-^ 5^ d,,d,, ■ ■ ■ di,{Zi, ■ ■ ■ Zi,r n - ^ifDn{f{z)){l + 0(1)). 

I<h<i2<—<ie<l+p l<j<k<e 

We now use Theorem 1.1 for Dn{f{z)). Noting, in particular, that G(l + z) = T{z)G{z) and 
P(0) = e^°, we obtain after a straightforward calculation that 

Dn{z'f{z)) = i-ir' Yl • • • ^^m- . . , Ai + 1, + 1, ■ ■ ■ , A, + 1, ■ ■ ■ ), 

l<il<i2<-<H<e+p 

where TZ is the r.h.s. of (1.12) where all (3j are replaced with Pj with the exception of Pj, j = 
ii, . . . which are replaced as indicated in the argument of TZ. Note once again that each sum is 
over indices in the range 1, ...,£ + p and we use a special numbering of indices (cf. (6.3)). Finally, 
recalling (6.2), we obtain 



(6.11) DM{z))={l[4- 



TZ{...,pi, + i,Pi, + i,...,pi, + i,..: 

I<ii<i2<---<ie<e+p 



which is the statement of Theorem 1.13 for V{z) analytic in a neighborhood of the unit circle. 
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6.2. Extension to smooth V{z). If V{z) is just sufficiently smooth, in particular C°°, on the unit 
circle C so that (1.14) holds for s from zero up to and including some s > 0, we can approximate 
V{z) by trigonometric polynomials V''"'\z) = Yl^~-p(n) ^kz^, z £ C. First, consider the case when 
maxj^fc \^f3j — ^Pk\ < 1- Then 2niaxj |3?/9j — uj\ < 1, where lo is defined by (4.63). We set 

(6.12) p=[n'--''], 1/ = 2max|3ft/3j -a;| +ei, 

j 

where ei > is chosen sufficiently small so that v <1 (square brackets denote the integer part). 

First, we need to extend the RH analysis of the previous sections to symbols which depend on 
n, namely to the case when F in / is replaced by V'^'^\ (We will denote such / by /(z, T^*^")), and 
the original one, by f{z,V).) We need to have a suitable estimate for the behaviour of the error 
term in asymptotics with n. For fixed /, our analysis depended, in particular, on the fact that 
f{z)~^z~"- is of order e"*^", £ > 0, for z G S""*^ (see Section 4.3), and similarly, f{z)~^z"' = 0{e~^'^) 
for z G S Here the contours S°"*, S are outside a fixed neighborhood of the unit circle 
(outside and inside C, respectively). If V is replaced by V^'^\ let us define the curve S outside 
^T=oUj by 

(6.13) ^=(^l + 7^je^^ 7>0, 



and E outside UJ^qUj by 

(6.14) ,= (i-^^I^) e 



p J 

Inside all Uj, the curves still go to zj as discussed in Section 4.2. Let the radius of all Uj be 
2'ylnp/p. We now fix the value of 7 as follows. Using the condition (1.14) we can write (here and 
below c stands for various positive constants independent of n) 



(6.15) 

/ P , „ . , 1/2 




^ i i;sT/i2l / (1 ± 37lnp/p)= 




\k=l 



(l±37lnfe/fc)±2M 1 



1/2 



< c 



k=l J \k=l 



1 + 



In^A;^ 



1/2 



where z G z £ dUj n {\z\ > 1} (with "+" sign in "±"), and z e z G dUj H {\z\ < 1} 

(with "— " sign). We now set 

(6.16) 37 = s-(l + £2)/2, £2>0, 
and then 

(6.17) |W")(z)|<c, |6+(z,"l/("))| < c, |6_(z,F('^))| < c, for ah n 

uniformly on E°"*, S ""^^j dUj's, and in fact in the whole annulus 1 — 37^ < jz] < 1 + 37^. 

It is easy to adapt the considerations of the previous sections to the present case, and we again 
obtain the expansion (4.60) for the jump matrix of R on dUj. Note that now |C(^)| = O(n'^lnn) 
and \z — Zj\ = In n/n^~'^ as n ^ 00 for z G dUj, and therefore using (4.55), (4.15), (4.17), (4.13), 
(4.9) and the definition of v in (6.12), we obtain, in particular, 

(6.18) n--3A,(^)n— 3 = q (;^^) , z G U]L,dUj. 
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Furthermore, as follows from (6.13), (6.14), (6.17), and (4.57,4.58), the jump matrix on and 
S is now the identity plus a function uniformly bounded in absolute value by 

,l-u \ 2maxj |3ff/3j| ^ In n ^ ^" 



(6.19) c 



Inn 



1 ± 7(1 - i^)^ ) < cexp {-|(1 - i/)n^ Inn} n'^i^'^'^^^j mi, 



where the upper sign corresponds to and the lower, to 

The RH problem for R{z) (sec Section 4.3) is therefore solvable, and we obtain R{z) as a series 
where the first term i?i is the same as before, and for the error term there holds the same estimate 
for z outside a fixed neighborhood of the unit circle, e.g., at 2; = 0. 

This, in particular, implies that the formulae (6.7,6.8) hold for f{z, F^")) (in /, we substitute /3j 
for Pj-. note that the condition < < 1 is satisfied). 

We will now show that replacing V^"^ with V in the symbol of the determinant Dn{f{z,V^^^)) 
results in a small error only, so that (6.10) still holds with V used in Dn{f{z, V)), and F^") in d/s 
and in D„(z^/(2;, F*^"))). Then, proceeding as before, we obtain the statement of the theorem for 
Dn{z^f{z,V(^))) as, by (1.14), 



(6.20) 



b±{zjM'''>) = b±{z„V) 



1 + 



1 



in (6.8). Recall a standard representation for a Toeplitz determinant with (any) symbol f{z): 

We have from this formula, (6.20), and Theorem 1.1 for Dn{\f{z, V)\) and Dn{f{z, V)), if s(l — i^) > 
1, 



(6.21) D,(/(z,y))-L>„(/(z,l^W) 



< 



——^ ... n \-''' - -'''\" n X ( 1 + -M^^-^^i - 1) 

^ > ■ -^0 \<n<k<n .7=1 



\<j<k<n J 



< C 



Er=o((si°.)'+w.)') 



< c 



(l-v)s-l 



Dnif{z,V)) 



n 



n 

-((i-^)^-i-Er=o((s>«,)2+m)')) 



Therefore, 
(6.22) 

DMizM""^)) = DM{z,V)) 1 1 + 



Dn{f{zM''^))-Dn{f{z,V)y 
Dnif{z,V)) 



Dr,{f{z,V)){l + o{l)), 



if 

(6.23) 



s > 



1 - u 



Under the condition (6.23) and the one under which Theorem 1.1 holds, e.g. C°° (see Remark 
1.7), we then obtain the statement of the theorem for Dn{z^f{z,V^"'^)) as mentioned above. The 
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theorem (with Remark 1.15) for Dn{z^ f{z,V)), and hence for Dn{f{z,V)), immediately follows 
from a similar to (6.21,6.22) analysis applied to 

n^(.'n., V)) ^ V""))) ( 1 - ^"'-'^''-„';':'»-°t'''"'^" ) ■ 

V D„{z'^f{z,VW)) J 

~ (r) 

The ratio in the brackets is o(l) under the condition (6.23) in which Pj are replaced by /3j (and 

the condition under which Theorem 1.1 holds). As ei can be arbitrary close to zero, this condition 
together with (6.23) (note that these conditions are consistent with (6.16) and the requirement that 
7 > 0) and (1.15) for Theorem 1.1 yield the estimate (1.28). 

7. Hankel determinants. Proof of Theorem 1.20 

Consider the Hankel determinant with symbol w{x) on [—1, 1] given by (1.36). In this section 
we will find its asymptotics using the relation to a Toeplitz determinant established in Theorem 
2.6. Let X = cos^, z = e*^, < ^ < tt. In particular, 

Xj = cos0j, Zj = e^^^, j = 0, 1, . . . , r + 1, = < < • • • < dr+i = tt. 

First, we find an even function / of the angle 9 related to w{x) by (1.39). The Toeplitz determinant 
D2n{f{z)) with this symbol enters the connection formula (2.20). Denote 

(7.1) 4 = e(2-%)^ i = 0,...,r + l. 
Then, recalling (1.6), note that 

(7.2) \x - Ajf = I cos ^ - cos ^jf 



e - 

2 sin sm ■ 



2-2«>|^_2^.|2".|2_2^.|2">, 



and 

(7.3) I sin6'| = 2~"^|z — ^olk — -2r+i|- 

We see that f{z) will have m + 1 = 2r + 2 singularities at the points zq = 1, Zr+i = —1, zj, Zj, 
j = l,...,r. 
Observe that 



j=i j=i j=i 



Note that f3o = Pr+i = and wc have the jumps with —(3j at zj and +f3j at z'j. In particular, 
the sum over all /3's is zero as noted in the introduction. Note that as 6j = 7r/2 — arcsin Aj, we have 
in (7.4) 

r I ^ 

(7.5) e''"" ^^'=1 zj^' z'f = exp 2i ^ arcsin A^- 
Collecting the above observations and denoting 

r+l 

(7.6) A = ^aj, 

j=0 

we have by (1.39) (where we single out a multiplicative constant for convenience) 

(7.7) f (z) = w{x)\ sin e\ = Cf{z), C = 2-^^-^ exp | 2i ^ /3j arcsin A^- j , 
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where 



Here V{eJ'^) = f/(cos6'). Thus f{z) is the symbol of type (1.2) with G (-1/2, 1/2]. Therefore, 
if ?R.Pj G (—1/2, 1/2), j = 1, . . . , r, we can apply Theorem 1.1 to D2n{f{z)), and obtain 

(7.9) L>2n(/(^)) =C2-L»2n(7(^)) = C2-exp (^nVo + f2kvi^ 6+(l)-^"°-i6+(-l)-^"'-+i-i 

f[ b+izj)-^^''^+^^h^{Zj)-^^''^-^^^ X (2n)2 2:;=lH-/3|)+(2«0+l/2)2+(2a.+i+l/2)2^^^) 



X 

i=i 

G(1 + g^. + /?^.)2g(i + g . - p^f G{1 + 2ao + 1/2)2 ^ 2a,+i + 1/2)^ 

y G'(l + 2a,)2 G(l + 4ao + l) G(l+4a,+i + l) ^^^>>^ 



3f?aj >-^, 3ft/3j- G f-^,^ ) , j = 0,l,...,r+l. 



2sin ^ 


-4:{ajak-l3j(3k) 


„ . 6*7 -1- 6*^; 

2sm — 


-i{ajak+(}jl3k) 


2 




2 





where we used the fact that by the symmetry of V{z), Vk = V-k, and hence b-^-{zj) = b-{zj), 
b-{zj) = b+{zj), 6+(±l) = 6_(±1). In the above expression for "P, m = 2r + 1, and the points are 

numbered as in Theorem 1.1, namely, zq = zq = 1, ao = 2ao + 1/2, Pq = /3o = 0; zj = zj, aj = aj, 
= -/^i' i = 1' • • • i'^; 2r+i = Zr+1 = -1, 3r+i = 2a,.+i + 1/2, J3r+i = f3r+i = 0; Zj = e^''''-z~^_^-^_j, 
otj ~ ^m+i-ji Pj = Pm+i-j, j = r + 2, . . . ,m. Expression for V can be written in terms of Xj. 
Namely, it is not difficult to obtain by induction that 

(7.10) V = 2~2(ao-(-ar-+i+l/4) JJ 

0<j<A;<r-+l 
r 

X JJ |2sin^j-p^H+^l+"^) X e2^(2^+l)E,"=i/3iarcsinA^g27riEo<,<fc<r-+i{"iA-"'=^.). 

Assume first that V{z) is analytic. To use (2.20), we need to calculate the asymptotics of the 
product 'l>2n(l)*I'2n(— 1)- In order to do this, consider Y^"'\z) as z — > zj in such a way that z G z{I), 
where z{I) is the pre-image in the z-plane of the / sector of the (^-plane (see Figure 2 and Section 
4.2). Tracing back the transformations of the RHP, we obtain 
(7.11) 

Yi-)(z)=nz) = Siz)(^^^^l_,^, fj={I + R^{\z))P,^iz)(^^^^l_,^^ J), zez{I), 

where the parametrix Pzj{z) at Zj is (see Section 4.2): 

(7.12) P,.{z) = E{z)^j{QFj{z)-''-^z'"'^/\ 
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with E{z) given by (4.47). Substituting all the expressions into (7.11), we obtain 

(7.13) Y(-\z) = {I + 4\z)Mzr (^^^ J) (c^^^.(-)-7"^T' 

^\ e^'^(^^+2"i)y'^^^^V^j(2)/(^)-^ Fj{z)J^ 

Note that the expansion of F~^{z) as z ^ zj is given by (4.51). Using that we further obtain 
for the last matrix in (7.13) 

(7.14) Fj{z)f{z)-^ = r/7ie^^(^^-°^^/2)zJ^u-"^(l + 0{u)), u = z-Zj, C e /. 
Thus, 

Fi(z)-^ \ ( e^™i 



To estimate ^(C) for C ^ (i-^-, ^ ^ -Zj), assume first that all ctj 7^ 0. Now substituting (4.34) 
into (4.32), dropping the terms of order u^"-' in the second column (we will denote thus modified 
Y{z) by Y^z)) we obtain the following limit for the combination needed in (7.13) (here tilde over 
the limit sign means that we have to drop tx^"-' terms before taking the limit): 

(7.16) lim„^o(^ Q gM/3,+2a,)J *i(C) g-i.a,J (e*""^^'^/^ "0"^ 

where 



(7.17) M 



r{a,+P,)r{l-aj-l3,)_ 

-iiraj 1 piTTOj r{l+aj+(3j) 

r(-f3j-aj) ^ r(aj-,(3j)r(l-aj+/3j) 

This expression can be simplified. Namely, the 11 matrix clement 

^ ^ r(/3,-a,)V' sin7r(/?,-a,); 

_ V{-2aj) sin(-27raj) _ r(l + - /3j) 
~ Vifij - aj) sin7r(/3j - aj) ~ r(l+2Q;j) 

r(i + «,+/?,) 



Similarly, 

(7.19) M21 = 



r(l + 2a,) • 



Thus 

(7.20) M 



/ r(l+aj-/3j) r(2«j) - 

/ r(l+2aj) r(a,+/3,) 
I r{l+aj+l3j) r{2aj) 



Substituting the just found limit and (4.53) for {V{z)l{C,^^F^{z)f into (7.13), we obtain 
(7.21) Y [z,)-{I+r. )L. , L- " (^_jv^^^^-i^-i„a,+/3, _Mi2M7V~"^+^^^7V ' 
where = R^\zj), and r/j, are given by (4.52,4.54). 



40 



P. DEIFT, A. ITS, AND I. KRASOVSKY 



(n) 

Note that the matrix Lj ' has the structure 



(7.22) 



where L depends on n only via the oscillatory terms 2". 
From (3.1) and (7.21) at = 1, 

$2n(l) = n?"^(l) = 43(1 + 0(n-2--'=/^'=-i)). 



(7.23) 



Prom (7.21,7.20,4.54,4.52), we obtain using the doubling formula for the F-function, 

F(l + a;) ^ 



(7.24) 



F(l + 2x) 22^F(x + l/2)' 



the following main term of $2n(l)- 
(7.25) L^^t = Mo,2i/xor?o-^(2n)2"o+V2 



0fg(Vb-^^{i))/2+iE-=i(T-e,)/3i ^ 



24ao+ir(l + 2ao) 
Similarly, we obtain 
(7.26) $2n(-l) = 4'+tn(l + 



2 sin — 
2 



-2a,- 



22(ao-ar+i)^2ao+l/2 



(2n) 
r+1,11 



24a.+i+lr(l+2ar+l) 



n 

5=1 



2cos^ 



-2a, 



2" 2(ao-ar+i)j^2ar+i-|-l/2 



Therefore 



(7.27) $2n(l)^2n(-l) = 



7re 



V^)-(V^(l)-fy(-l))/2-|-2j X;^=i arcsin A,- 



24(ao+a.+i)+2r(i + 2ao)r(l + 2a^+i) 



X JJ|2sinej|-2°%2(ao+a,+i)+l^;L^(5^^-2max,A-l^^_ 

Substituting (7.9) and (7.27) into (2.20) we obtain (1.37) squared. We use the following obser- 
vations in the process: 

• Since T4 = V-k, 

6+(±l) = e(^(±i)-^°)/^ 

• The following elementary identity holds 



(7.28) n 

0<j<k<r+l 



r, . - (^k 

2sm — 


-{otjotk-pjpk) 


9j + 9k 
2sm — 


-{ajak+PjPk) 


2 




2 





2-^0<j<k<r+l°'j°'k JJ |_\^. _ ;^^|-("jQ;fe+/3j/3fc) 

0<i<fc<r+l 



A,A,-1 + ^(1-Aj)(l-A|) 



(7.29) 



• Applying the doubling formula (2.39) we easily obtain that 

2 - G(l/2)2 



G(l + 2a + 1/2)2 _ 
G'(l + 4a + l) ^^^ + ^^)-^ ^ G(l + 2a)2- 
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If V{z) = Vr{z) is real-valued for z G C, and a.j e M, i(3j € M, j = 0, ...,m, then the 
weight /(-z), z G C, is positive, and therefore Dn{w) is positive. Then (1.37) represents the 
correct branch of the square root. Since Dn{w) is continuous in ctj, (3j, and the parameter h in 
F(2;) = Vr{z) + {y{z) — Vr{z))h, h G [0, 1], and the error term is uniform in these parameters (see 
Section 4.3), the formula (1.37) has the correct sign in general. This finishes the proof for analytic 
V{z). The extension to smooth V{z) is carried out similarly to the argument in the previous section 
by using the standard multiple- integral representation of a Hankel determinant. 
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